A SHAPE THEOREM FOR THE SPREAD OF AN INFECTION 



Abstract. In [KSb] we studied the following model for the spread of a rumor or 
infection: There is a "gas" of so-called A-particles, each of which performs a con- 
tinuous time simple random walk on Z d , with jumprate D^. We assume that "just 
before the start" the number of A-particles at x, Na(x, 0— ), has a mean /x^ Poisson 
distribution and that the Na(x,0—), x £ Z d , are independent. In addition, there are 
_B-particles which perform continuous time simple random walks with jumprate Dg. 
We start with a finite number of B-particles in the system at time 0. The positions 
of these initial B-particles are arbitrary, but they are non-random. The B-particles 
move independently of each other. The only interaction is that when a B-particle and 
an A-particle coincide, the latter instantaneously turns into a B-particle. [KSb] gave 
some basic estimates for the growth of the set B(t) := {x £ Z d : a B-particle visits x 
during [0,*]}. In this article we show that if D A = D B , then B(t) = B(t) + [-|, \] d 
grows linearly in time with an asymptotic shape, i.e., there exists a non-random set 
Bq such that (l/t)B(t) — > Bq, in a sense which will be made precise. 



1. Introduction. 

We study the model described in the abstract. One interpretation of this model 
is that the S-particles represent individuals who are infected, and the A-particles 
represent susceptible individuals; see [KSb] for another interpretation. Bit) repre- 
sents the collection of sites which have been visited by a S-particle during [0, i], 
and B(t) is a slightly fattened up version of Bit), obtained by adding a unit cube 
around each point of Bit). This fattened up version is introduced merely to sim- 
plify the statement of our main result. It is simpler to speak of the shape of the 
set (l/t)B(t) as a subset of M d , than of the discrete set (l/t)B(t). 

The aim of this paper is to describe how the infection spreads throughout space as 
time goes on. In [KSb] we proved a first result in this direction in the case Da = Dg. 
We proved that under this condition there exist constants < C<i < C\ < oo such 
that almost surely 
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where 

C(r):=[-r,r] d . (1.2) 

(1.1) gives upper and lower bounds which are linear in time, for B(t), the region 
which has been visited by the infection during [0, t]. However, the upper and lower 
bounds in (1.1) are not the same. The principal result of this paper is a so-called 
shape theorem which gives the first order asymptotic behavior of the region B(t). 
It shows that (l/t)B(t) converges to a fixed set B . Thus, not only is the growth 
linear in time, but B(t) looks asymptotically like (a scaled version of) Bq. This of 
course sharpens (1.1) by 'bringing the upper and lower bound together'. However, 
the result (1.1) is a crucial tool for proving the shape theorem. We do not know 
of a shortcut which proves the shape theorem without much of the development of 
[KSb] for (1.1). The precise form of the shape theorem here is as follows: 

Theorem 1. Consider the model described in the abstract. If Da = Db, then 
there exists a non-random, compact, convex set Bq such that for all e > almost 
surely 

(1 - e)B C -B(t) C (1 + e)B for all large t. (1.3) 

t 

The origin is an interior point of B , and B is invariant under reflections in 
coordinate hyperplanes and under permutations of the coordinates. 

Remark 1. It follows immediately from Theorem 1 and Proposition B below that 
the particle distribution at a large time t looks as follows: The numbers of particles, 
irrespective of type, that is Na(x, t) + Nb(x, t), x G Z d , is a collection of i.i.d. mean 
[iA Poisson variables plus a finite number of particles which started at time zero 
at fixed locations (these are the particles added as B-particles at the start). For 
every s > there are almost surely no A particles in (1 — e)tB and no S-particles 
outside (1 + e)tB for all large t. 

Shape theorems have a fairly long history and have become the first goal of many 
investigations of stochastic growth models. To the best of our knowledge Eden (see 
[E]) was the first one to ask for a shape theorem for his celebrated 'Eden model'. 
The problem turned out to be a stubborn one. The first real progress was due to 
Richardson, who proved in [Ri] a shape theorem not only for the Eden model, but 
also for a more general class of models, now called Richardson models. In these 
models one typically thinks of the sites of Z d as cells which can be of two types (for 
instance B and A or infected and susceptible). Cells can change their type to the 
type of one of their neighbors according to appropriate rules. One starts with all 
cells off the origin type A and cell of type B at the origin and tries to prove a shape 
theorem for the set of cells of type B at a large time. An important example of such 
a model is 'first-passage percolation', which was introduced in [HW] (this includes 
the Eden model, up to a time change). A quite good shape theorem for first-passage 
percolation is known (see [Ki], [CD], [Ke]). In more recent first-passage percolation 
papers even sharper information has been obtained which gives estimates on the 
rate at which (l/t)B(t) converges to its limit Bq (see [Ho] for a survey of such 
results). 
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Shape theroems for quite a few variations of Richardson's model and first-passage 
percolation have been proven (see for instance [BG] and [GM]), but as far as we 
know these are all for models in which the cells do not move over time, with one 
exception. This exception is the so-called frog model which follows the rules given 
in our abstract, but which has Da = 0, i.e., the susceptibles or type A cells stand 
still (see [AMP] and [RS] for this model). The present paper may be the first one 
which allows both tyes of particles to move. 

In nearly all cases shape theorems are proven by means of Kingman's subadditive 
ergodic theorem (see [Ki]). This is also what is used for the frog model. For this 
model one can show that the family of random variables {T xy } is subadditive, were 
T XyV is a version of the first time a particle at y is infected, if one starts with one 
infected particle at x and one susceptible at each other site. More precisely, the 
T XjV can all be defined on one probability space such that T XjZ < T XjV + T y 2 for all 
x,y,z G Z d and such that their joint distribution is invariant under translations. 
Unfortunately this subadditivity property is no longer valid if one allows both types 
of particles to move. Nevertheless, subadditivity methods are still heavily used in 
the proof of Theorem 1. However, we now use subadditivity only for certain 'half- 
space' processes which approximate the true process. Moreover, these half-space 
processes have only approximate superconvolutive properties (in the terminology 
of [Ha]). There is no obvious family of random variables with properties like those 
of the T Xj y. One only has some relation between the distribution functions of the 
H (t, u) for a fixed unit vector u, where H(t, u) is basically the maximum of (x, u) 
over all x which have been reached by a S-particle by time t ((x,u) is the inner 
product of x and u; for technical reasons H (t, u) will be calculated in a process 
in which the starting conditions are slightly different from our original process). 
These properties are strong enough to show that for each unit vector u there exists 
a constant \{u) such that almost surely 



Thus the S-particles reach in time t half-spaces in a fixed direction u at distances 
which grow linearly in t. Except in dimension 1, it then still requires a considerable 
amount of technical work to go from this result about the linear growth of the 
distances of reached half-spaces to the full asymptotic shape result. We will give 
more heuristics before some of our lemmas. 

Remark 2. Our proof in [KSb] shows that the right hand inclusion in (1.1) remains 
valid for arbitrary jumprates of the A and the S-particles. However, it is still not 
known whether the left hand inclusion holds in general. The lower bound for B(t) is 
known only when Da = Db, or when Da = 0, that is, when the A and S-particles 
move according to the same random walk (see [KSb]), or in the frog model, when 
the A-particles stand still (see [AMP],[RS]). 

Here is some general notation which will be used throughout the paper. \\x\\ 
without subscript denotes the £°°-norm of a vector x = (x(l), . . . , x(d)) G M d , i.e., 



lim -H(t, u) = X(u) 

n^oo t 



(1.4) 



xll = max \x(i 

Ki<d 
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We will also use the Euclidean norm of x; this will be denoted by the usual ||x||2- 
(x, it) denotes the (Euclidean) inner product of two vectors i,«6 M d , and denotes 
the origin (in Z d or M, d ). For an event £, £ c denotes its complement. 

Ki,K2,... will denote various strictly positive, finite constants whose precise 
value is of no importance to us. The same symbol Ki may have different values 
in different formulae. Further, Cj denotes a strictly positive constant whose value 
remains the same throughout this paper, a.s. is an abbreviation of almost surely. 
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2. Results from [KSb]. 

Throughout the rest of this paper we assume that 

D A = D B (2.1) 

and we abbreviate their common value to D. We begin this section with some 
further facts about the setup. We concentrate on the special case Da = D b . More 
details can be found in Section 2 of [KSb] which deals with the construction of 
our particle system. {St}t>o will be a continuous time simple random walk on Z d 
with jumprate D and starting at 0. To each initial particle p is assigned a path 
{71^4 (t, p)}t>o which is distributed like {S t }t>o- The paths tta(-,p) for different p's 
are independent and they are all independent of the initial Na(x, 0—), x G Z d . The 
position of p at time t equals 7r(0, p) +nA(t, p), and this can be assigned to p without 
knowing the paths of any of the other particles. The type of p at time s is denoted 
by 77 (s, p). This equals A for < s < 6(p) and equals B for s > 9(p), where 9(p), the 
so-called switching time of p, is the first time at which p coincides with an initial B- 
particle. Note that this is simpler than in the construction of [KSb] for the general 
case which may have Da ^ Db- In that case we had simple random walks {S v }t>o 
with jumprate D v for 77 G {A, B}, and there were two paths associated with each 
initial particle p : tt v (-, p), rj G {A, B}, with {7T v (t, p)} having the same distribution 
as {S^}. If p had initial position z, its position was then equal to z + ^(0, p) until 
p first coincided with a B-particle at time 9{p); for t > 9(p) the position of p was 
z + ita(9(p), p) + [ns(t, p) — ■kb{9{p)- i p)]. This depends on 9{p) and therefore on the 
movement of all the other particles. In the present case we can take tyb = ^ 'a, which 
has the great advantage that the path of p does not depend on the paths of the other 
particles. This is the reason why the case Da = Db is special. We proved in [KSb] 
that on a certain state space E , the collection of positions and types of all particles 



SHAPE THEOREM FOR SPREAD OF AN INFECTION 



5 



at time t, with t running from to oo, is well defined and forms a strong Markov 
process with respect to the cr-fields Tt = f^hyo^t+h^ * — ^> wnere is the cr-field 
generated by the positions and types of all particles during [0,£]. The elements of 
these (T-fields are subsets of E[°'°°), where E = n fc >i ((Z d U d k ) x {A, B}) . £[°'°°) is 
the pathspace for the positions and types of all particles. More explicit definitions 
are given in [KSb] but are probably not needed for this paper. It was also shown 
in [KSb] that if one chooses the number of initial A-particles at z, with z varying 
over Z d , as i.i.d. mean ha Poisson variables, then the process starts off in E and 
stays in Eo forever, almost surely. 

We write N v (z 7 t) for the number of particles of type r\ at the space-time point 
(z,t), z E Z d ,i] E {A,B}, while Na(z, 0— ) denotes the number of A-particles to 
be put at z "just before" the system starts evolving. Note that our model always 
has only particles of one type at each given site, because an A-particle which meets 
a B-particle changes instantaneously to a S-particle. Thus, if Na(z, 0— ) = N for 
some site z and we add M > S-particles at z at time 0, then we have to say that 
N A (z, 0) = 0, N B {z, 0)=N + M. 

We shall rely heavily on basic upper and lower bounds for the growth of B(t) 
which come from Theorems 1 and 2 in [KSb]. 

Theorem A. If Da = Db, then there exist constants < C2 < C\ < 00 such that 
for every fixed K 

P{C(C 2 t) c B(t) C C(2d*)} > 1 - ^ (2.2) 

for all sufficiently large t. 

We also have some information about the presence of A-particles in the regions 
which have already been visited by S-particles. The following is Proposition 3 of 
[KSb]. 

Proposition B. If Da = Db, then for all K there exists a constant C3 = Cs(K) 
such that 

P {there is a vertex z and an A-particle at the space-time point (z,t) while 
there also was a B -particle at z at some time < t-C 3 [tlogt] 1/2 } 

< for all sufficiently large t. (2.3) 

Consequently, for large t 

P{at time t there is a site in C (6*2^/2) which 

2 

is occupied by an A-particle} < (2.4) 

Finally we reproduce here Lemma 15 of [KSb] which gives an important mono- 
tonicity property. We repeat that in the present setup, with the Na{x,0— ) i.i.d. 
Poisson variables, our process a.s. has values in Eo at all times (see Proposition 5 
of [KSb]). 
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Lemma C. Assume Da = Db and let G So- Assume further that lies 
below o~( 2 ) in the following sense: 

for any site z G 7L d , all particles present in at z are also present in at z, 

(2.5) 

and 

at any site z at which the particles in a*- 2 -* have type A, 

the particles also have type A in (2-6) 

Let 7Ta(-, p) be the random walk paths associated to the various particles and assume 
that the Markov processes {Y^} and {Y^} are constructed by means of the same 
set of paths 7r^(-, p) and starting with state a^ 1 ' and , respectively (as defined in 
Section 2 of [KSb], but with tta(-,p) = ^b(',p) f or all s,p; see (2.6), (2.7) there). 
Then, almost surely, {Y^} and {y} 2 ^ 1 } satisfy (2.5) and (2.6) for all t with 
replaced by Y t , i = 1,2. In particular, G E . 

In particular, this monotonicity property says that if is obtained from > 
by removal of some particles and/or changing some B-particles to A-particles, then 
the process starting from has no more S-particles at each space-time point 
than the process starting from o~( 2 K We note that this monotonicity property holds 
only under our basic assumption that Da = Db- 

3. A subadditivity relation. 

In this section we shall prove a basic subadditivity relation and deduce from 
it that the S-particles spread in each fixed direction over a distance which grows 
asymptotically linearly with time. This statement is ambiguous because we haven't 
made precise what "spread in a fixed direction" means. Here this will be measured 
by 

max{(i, it) : x G B(t)}, (3.1) 

where u is a given unit vector (in the Euclidean norm) in W 1 (see the abstract for 
B). In addition we will not prove subadditivity (which is an almost sure relation), 
but only superconvolutivity, in the terminology of [Ha] (which is a relation between 
distribution functions). The tool of superconvolutivity in other models with no 
obvious subadditivity in the strict sense goes back to [R], and was also used in [BG] 
and [W]. 

We define the closed half-space 

S(u,c) = {x G R d : (x,u) > c}. 

Given a u G S d ~ x and r > we consider the half -space process corresponding to 
(u,—r) (also called (u,—r) half-space-process). We define this to be the process 
whose initial state is of the form 

N A (x, 0-) = if x S(u, -r) and the N A (x, 0-), x G S(u, -r), 

are i.i.d., mean pa Poisson variables. 
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In addition the particles at xo : - r are turned into B-particles at time 0, where £o,-r 
is the site in S(u, —r) nearest to the origin (in £°°-norm) with Na{xq- t , 0—) > 0; 
if there are several possible choices for xo t - r , the tie is broken according to some 
deterministic rule chosen in advance. There will be many other occasions were 
ties may occur. These will be broken in the same way as here, but we shall not 
mention ties or the breaking of them anymore. Note that no extra S-particles are 
introduced at time 0, but that only the type of the particles at xq- t is changed. 
Thus, 

N A (x, 0) + N B (x, 0) = N A (x, 0-) for all x. (3.2) 

From time on the particles move and change type as described in the abstract. 
Note that only the initial state is restricted to S(u, — r). Once the particles start 
to move they are free to leave S(u, — r). The (u, — r) half-space process will often 
be denoted by V (u, —r). 

We further define the (u, — r) half-space process starting at (x,t). This process 
is defined for times t' > t only. We define it as follows: at time t let x 0j - r (t) be the 
nearest site to x which is occupied in the (u, — r) half-space process. We then reset 
the types of the particles at xo- r (i) to B and the types of all other particles present 
in the (u, — r) half-space process at time t to A. The particles then move along the 
same path in the (w, — r) half-space process starting at (x, t) as in V h (u, — r) (which 
starts at (0, 0)). However, the types of the particles in the (w, — r) half-space process 
starting at (x,t) are determined on the basis of the reset types at time t. Thus 
the half-space process starting at (x,t) has at any time only particles which were 
in iS(w, — r) at time 0. Moreover, at any site y and time t' > t, V (u, — r) and 
the (w, — r) half-space process started at (x,t) contain exactly the same particles. 
We see from this that the paths of the particles in the (-u, — r) half-space processes 
starting at (x,t) and at (0,0) are coupled so that they coincide from time t on, 
but the types of a particle in these two processes may differ. Lemma C shows 
that if there is a S-particle in V (u, — r) at x at time t, then in this coupling any 
B-particle in the (w, — r) half-space process starting at (x,t) also has to have type 
B mV h {u,-r). 

The coupling between the two half-space processes clearly relies heavily on the 
assumption Da = Db, so that we can assign the same path to a particle in the 
two processes, even though the types of the particle in the two processes may be 
different. 

It is somewhat unnatural to start the (u, — r) half-space process with S-particles 
at xo,- r in case r < 0, so that the origin does not lie in the half-space S(u, —r). 
We shall avoid that situation. We can, however, use the (u, — r) half-space process 
starting at (x,t). This is well defined for all r. We merely need to find the site 
nearest to x which has at time t a particle which started in S(u, — r) at time 0. We 
can then reset the type of the particles at this site to B at time t. We shall consider 
the (w, — r) half-space process starting at (x, t) mostly in cases where we already 
know that x itself is occupied at time t in the (u, — r) half-space process. 

Finally we shall occasionally talk about the full-space process and the full-space 
process starting at (x,t). These are defined just as the half-space processes, but 
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with r = oo. In particular, the full-space process starts with S-particles only at the 
nearest occupied site to the origin and (3.2) applies. The full-space process starting 
at (x,t) has S-particles at time t only at the nearest occupied site to x. The type 
of all particles at other sites are reset to A at time t. By stationarity in time, the 
full-space process started at (x, t) has the same distribution at the space-time point 
(x + y,t + s) as the full-space process (started at (0, 0)) at the point (y, s). Again 
we shall use the same random walk paths tva for all the full state processes and the 
half-space processes, so that these processes are automatically coupled. We shall 
denote the full-space process by . 

We point out that if < T\ < r 2 , and if ||xo,-r|| < fi/Vd, then £o,-r £ 
S(u, — ri) C S(u, — r 2 ). In this case, both V h (u,—ri) and V h {u,—r2) start with 
changing the type to B at the site xo t - r only. By Lemma C, at any time 

any S-particle in V h (u, — r{) is also a B-particle in V h (u, — r 2 ). (3-3) 

This comment also applies if V h {u, — r 2 ) is replaced by V* (which is the case r 2 = 
oo). 

It seems worthwhile to discuss more explicitly the relation of the full-space pro- 
cess to our process as described in the abstract. The latter has some S-particles 
introduced at time at one or more sites, in addition to the Poisson numbers of 
particles, Na(x, 0— ), x G Z d . If exactly one S-particle is added at time 0, and this 
particle is placed at 0, then we shall call the resulting process the original process. 

Suppose we want to estimate P{A(xo)} in the full-space process, where 

xq '■= the nearest occupied site to the origin at time in , (3.4) 

A is some event and A(x) is the translation by x of this event (which takes Na(0 7 s) 
to N A (x,s)). Then, for C a subset of Z d , 

P{x eC,A(x ) in V f }= J2 p i 

x = x,A(x)} 

xec 

< P{x is occupied at time 0, A(x) in V^} 



xec 



oo 



= ^2^2e~^ A ^-P{A\there are k S-particles at at time 0}. 

xeCk=i k - (3.5) 

(The probability in the last sum is the same in 

V f 

as in the original process.) On 

the other hand, in the original process we have 
P{A in original process} 

= X^ e ~ ffc p-j P{.4j there are k S-particles at at time 0}. 

fc=i ^ '' (3.6) 
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Comparison of the right hand sides in (3.5) and (3.6) yields the crude bound 

P{xq G C, A(xq) in the full-space process} 

< (cardinality of C)/j,aP{A in original process}. (3.7) 

We shall repeatedly use a somewhat more general version of this inequality (see for 
instance (3.25), (3.77), (3.78), (5.33)). Suppose s > is fixed and X is a random 
vertex in Z d , and suppose further that 

P{in V f , A(X) but (X, s) is not occupied} = 0. (3.8) 

(Note that this is satisfied if (X, s) is occupied almost surely in V* .) Let C C 7L d as 
before. Now, given that there are k > 1 particles at the (non-random) space-time 
point (x, s), the full-space process starting at (x, s) is simply a translation by (x, s) 
in space-time of the original process, conditioned to start with k — 1 points at the 
origin and one B-particle added at the origin. Therefore, essentially for the same 
reasons as for (3.7), 

P{X G C, A(X) in the full-space process starting at (X, s)} 

< (cardinality of C)haP{A in original process}. (3.9) 

For a rather trivial comparison in the other direction we note that if P{A mV^} = 
for the full-space process, then we certainly have for each k > 1 that 

= P{A in , xq = 0, k particles at xq} 

= P{A mV f ,k particles at 0} 

\ UA ] k 

= e" MA ^-^P{^l|there are k S-particles at at time 0}. (3.10) 
k\ 

This implies, via (3.6) that also P{A in V f } = 0. 

It is somewhat more complicated to compare with the process described in 
the abstract if more than one S-particle is introduced at time 0. Rather than 
develop general results in this direction we merely show in our first lemma that it 
suffices to prove (1.3) for the full-space process. 

Lemma 1. // (1-3) holds in , then it also holds in the original process of the 
abstract with any fixed finite number of B -particles added at time 0. 

Proof. The preceding discussion shows that if (1.3) has probability 1 in V* , then 
it has probability 1 in the original process (with one particle added at the origin at 
time 0). By translation invariance (1.3) will then have probability 1 in the process 
of the abstract with one particle added at any fixed site at time 0. 

Lemma C implies that one can couple two processes as in the abstract, with 
collections of S-particles C A^ added at time 0, respectively, in such a way 
that the process corresponding to A^ always has no more S-particles than the 



10 HARRY KESTEN AND VLADAS SIDORAVICIUS 

one corresponding to A^ 2 \ Therefore, if the left hand inclusion in (1.3) holds when 
only one B-particle is added at time 0, then it certainly holds if more than one 
B-particle are added. 

It follows that we only have to prove the right hand inclusion in (1.3) for the 
process from the abstract with more than one particle added, if we already know it 
when exactly one particle is added. Assume first that we run this last process with 
one .B-particle po added at zq. We now have to refer the reader to the genealogical 
paths introduced in the proof of Proposition 5 of [KSb] . The right hand inclusion 
in (1.3) then says that for all e > 

P{there exist genealogical paths from zq to some point 

outside (1 + e)tB for arbitrarly large t} = 0. (3.11) 

From the construction of the genealogical paths in Proposition 5 of [KSb] and the 
fact that a.s. there are only finitely many S-particles at finite times (see (2.18) in 
[KSb]) it is not hard to deduce that 

{B{t) <f_ (1 + e)tBo at time t if one adds a S-particle pi 
at Zi, 1 < i < k, at time 0} 

= {there is a genealogical path from some z^, 1 < i < k, 
to the complement of (1 + s)tBo at time t if one 
adds a B-particle pi at z^, 1 < i < k, at time 0} 

k 

C {there is a genealogical path from Zi to the complement of 

i=i 

(1 + e)tB at time t if one adds a S-particle pi at zi at time 0} 

(3.12) 

(the Zi do not have to be distinct here). It follows that 

P{B(t) <f_ (1 + e)tB for arbitrarily large times t if one 

adds a S-particle pi at Zi, 1 < i < k, at time 0} 

k 

< ^P{there are genealogical paths from Z{ to the complement of (1 + e)tB 
i=i 

at arbitrarily large times t if one adds a S-particle pi at Zi at time 0} 
= (by (3.11)). 

Thus the right hand inclusion in (1.3) holds a.s., even if one adds k S-particles at 
time 0. ■ 

We recall that 

V (it, — r) is short for the (u, — r) half-space process, 
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V* is short for the full-space process, 
and we further introduce 

B h (y, s; u, —r) := {there is a B-particle at (y, s) in V h (u, —r)}, 

h(t,u,—r) = max{(j,«) : B h (x,t;u, — r) occurs}. 



(3.13) 
(3.14) 



P or will denote the probability measure for the original process (with one S-particle 
added at the origin at time 0); E or is expectation with respect to P or . (The 
superscripts h, f and or are added to various symbols which refer to a half-space 
process, the full-space process, or the original process, respectively). We use P 
without superscript if it is clear from the context with which process we are dealing 
or when we are discussing the probability of an event which is described entirely in 
terms of the Na(x, 0— ) and the paths tta- 

The following technical lemma will be useful. It tells us that, with high proba- 
bility, V h (u, — r) moves out in the direction of u at least at the speed C4, provided 
r is large enough (see (3.15)). Its proof would be nicer if we made use of the fact 
that even the (it, 0) halfspace-process has, with a probability at least 1 — t~ K , B- 
particles at time t at sites x with (x,u) > Ct, for some constant C > 0. However, 
it takes some work to prove this fact and we decided to do without it. The lemma 
itself is proven by recursively constructing a sequence of space-time points which 
move out in the direction of u along an exponentialy growing sequence, so that 
there is only an exponentially small (in k) probability that the k-th point is not 
occupied in the (u, —r) half-space process. 



Lemma 2. Let C±, C2 be as in Theorem A and let 



2\fdC 1 C 2 
32Vdd + C 2 ' 



(3.15) 



For all constants K > 0, there exists a constant ro = ro(K) > such that for 
r > r 



P< h(t, u, — r) < CaX for some t >t\ := 



1 



1 + 



Co 



Asfdd L 32^^ 



r} < r~ K . (3.16) 



Proof. 

Step 1. For k > 1 define the times 



tk 



1 



1 + 



Co 



and the real numbers 



dk 



Co 



1 + 



Co 



32 v / dCi L 32Vdd J 
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Also define for each k > 1 the event 

T>k '■= {B h (xk, tk] u, — r) occurs for some Xk which 

satisfies (xk,u) > dk and ||xjfc|| < 2C\tk\- (3-17) 

In this step we shall reduce the lemma to an estimate for the probability that U k 
fails for some k > 1. Indeed, assume that V k occurs for all k > 1. By definition, 
there is then a B-particle at (xk,tk) in the (u, —r) half-space process (starting at 
(0,0)), so that 



Co 



Co 



k 



h(t k ,u, -r) > <*»,«> > 4 = + r, *> 1. (3.18) 

Recall that jF t is defined in the beginning of Section 2. In addition to (3.18), we 
have on the event {{xk, u) > dk}, for k > 1, 

P{h(t,u,-r) < ^d k for some t G tfc+i)|^t fe } 
< P{each S-particle in V h (u, — r) at (xk,tk) moves during 

[tk,tk+i] to some site x with (x,u) < -dk} 

<P{ min (S q ,u) < --d k = -C 4 t k+ i} 
q<tk+i—tk I 

<K ie xp[-K 2 t fc+1 ] (3.19) 

for some constants Ki, K 2 depending on d, Da only; see (2.42) in [KSa] for the last 
inequality. It follows that the left hand side of (3.16) is bounded by 

oo 

P{V k fails for some k > 1} + K x exp[-K 2 t k ]. (3.20) 

k=i 

Step 2. We shall now derive a recursive bound for f^i<j<k^j- Assume that 
\~\i<j<k-\Dj occurs for some k > 2. Consider now the full-space process starting 
at (xk-i, tk-i)- Define the following events for this process: 

£fc,i := { & t time tk all occupied sites in 

Xk-i + C((C2/2)(tfc — tk-i)) contain in fact a S-particle}, 
£k,2 '■= {at time tk there is an occupied site in 

Xk-i + (C 2 /4)(t k - t k -i)u + C([logt fc ] 2 )}, 
£ k ,3 ■= {all particles in x k -i +C(2Ci(t k - t k -i)) 

at time tk-i started at time in S(u, — r)}, 
£k,4 '■= {there is no S-particle outside Xk-i + C{C\{tk — tk-i)) 

during [t k -i,t k ]}, 
£k,5 '■= {no particle which is outside Xk-i + C(2C\(t k — tk-if) 

at time t k -i enters x k -i +C(Ci(t k -t k -i)) during [t k -i, tk}}, 
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We claim that on 

Pfc-iH p| £ k ,i (3.21) 

l<i<5 

also V k occurs, provided r > some suitable ?r, independent of k, and k > 2. We 
merely need to make sure that v^[logtfc] 2 < (C2/8) (t k — tk-i) whenever r > r\. 
To prove our claim when k > 2, observe first that the occurrence of £k,i H £^2 
guarantees that at time tk there is a B-particle at some x k in Xk-i + (C f 2/4)(t i t — 
tk-i)u + C([logtfc] 2 ) C + C((C 2 /2)(t fc —tk-i))- Such a particle automatically 
satisfies 

C2 C2 

(xfc,u) > (x k -i,u) + — (t k -t k - 1 )-\Q[logt k ] 2 > d k -i + — {t k -tk-i) =d k . (3.22) 

It also satisfies \\xk\\ < 2C\t k , because ||xjfc_i|| < 2Cit k -i, and on £k,2, \\x k \\ < 
||xjfe_i|| + (C 2 /4)(£ fc - tfc-i) + [logtfc] 2 , while C 2 < C\. This particle at (x k ,t k ) 
is a B-particle in the full-space process starting at (xfc_i, tfe-i)- We are going to 
show that, in fact, it is also a S-particle in the (u, — r) half-space process starting 
at (x k -i, tk-i)- This will prove our claim, because the monotonicity property of 
Lemma C implies that any S-particle in the (u, — r) half-space process starting 
at (xk-i, tk-i) is also a S-particle in the (u,—r) half-space process (starting at 
(0, 0)), provided that there is a S-particle at (x k -i, tk-i) in the (u, — r) half-space 
process. (Note that this proviso is satisfied by the induction assumption that T> k -i 
occured.) We first observe that the particle at (x k , t k ) must at time t k -i have been 
in x k -i +C(2C 1 (t k - t k -i)), because x k e x k -i +C((C 2 /2)(t k - t k -i)) C x k - X + 
C{[C\/2){t k — tk-i)) and £k,5 occurs. By virtue of £^3 this particle then belongs 
to V h {u,—r) as well as to the (u,—r) half-space process starting at (xfc_i, tfe-i)- 
We still have to show that this particle also has type B in the (it, — r) half-space 
process starting at (x k -i, tk-i)- To this end we note that the particles starting 
outside Xk-i + C(2Ci(t k — t k -i)) at time tk-i do not influence the type of any 
particle at time t k in the full-space process starting at (x k -i, t k -i). Indeed, in this 
process the particles outside x k -i+C(2Ci(t k — t k _i)) start as A-particles, and since 
£fc,4 H £ k ,5 occurs, these particles do not meet any S-particle at or before time t k . 
Thus, whether the particle at (x k ,t k ) is also a S-particle in the (u, — r) half-space 
process starting at (xk-i,tk-i) depends only on the paths of the particles which 
were in x k -\ +C(2Ci(t k — tk-i)) at time t k -i (compare the lines following (2.37) in 
[KSb]). All these particles were particles in V h (u, — r) at time t k -i (on £ k ^), and 
hence also are in this half-space process at time t k . Thus the type of the particle 
at (x k ,t k ) is the same in the full-space process starting at (x k -i, t k -i) and in the 
(u, — r) half-space process starting at (x k -i, t k -i)- This justifies our claim that V k 
occurs for k > 2. We leave it to the reader to make some simple modifications in 
the above argument to show that T>i occurs on 

l<i<5 

where 

t = and V = {\\x \\ < if 3 log r}, 
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provided r\ is chosen large enough; xq is defined in (3.4) and K 3 is chosen right 
after (3.25) and depends on K, d and \ia only. 

We have now shown that on the event (3.21) also T>k occurs. If this is the case 
and also f\\<i<<s£k+i,i occurs, then T>k+i occurs etc. Consequently, for r > ri, 

P{V occurs, but some V k fails} 

5 

some xq with ||a;o|| < K3 logr, xo is occupied but £\^ fails} 

i=l 

00 5 

+ ^^^^P{for some x k -i with < 2C\tk-\ and (xk-i,u) > d^-i, 

k=2 i=l 

B h (xk-i,tk-i',u, — r) occurs, but £ k ,i fails}. (3.23) 

Step 3. In this step we shall give most of the estimates for the terms in the right 
hand side here for k > 2. The basic inequalities remain valid for k = 1 by trivial 
modifications which we again leave to the reader. For the various estimates we have 
to take all tk large. This will automatically be the case if r is large; we shall not 
explicitly mention this in the estimates below. 

We start with the estimate for the failure of £k,i- If £k,i fails, for a given 
(xk-i, tfc-i), then there must be some y G Xk-i +C((C2/2)(tk —tk-i) such that y is 
occupied by an A-particle at time tk in the full-space process started at (xk-i, tk-i)- 
Recall that if we shift the full-space process starting at (x,t) by (—x, —t) in space- 
time, then we obtain the full state process starting at (0,0). Moreover, if we 
condition on the event that x is occupied at time t, then, after the shift by (—x, —t) 
the NA(y, 0), y 7^ 0, are i.i.d. mean [ia Poisson random variables. Therefore, by 
summing over the possible values for Xk-i, 

P{for some Xk-i with ||xfc_i|| < 2C\tk-i and (xk-i,u) > dk-i, 
B h (xk-i, tfc-i; u, — r) occurs, but £k,i fails} 

occurs and in the full-space process started 

||x||<2Cit fc _i 

at (x,tk-i) there is an A-particle in x + C((C 2 /2)(t fc — t k -i)) at time t k } 
< P{0 is occupied at time and in there is an A-particle 

||*||<2Cit fc _i 

in C((C 2 /2)(tk - t k -i)) at time t k - t k -i}. 

To the right hand side here we can apply (3.7) (with C = {0}). This shows that 
the right hand side is at most 

K±\tk-i] d ^aP° t ^{at time tk—tk-i, there is an A-particle in C ((C2 /2)(tk — tk-i))} ■ 

The probability in the right hand side here is calculated for the original process 
with one particle added at at time 0. By (2.4) (with K replaced by K + d + 2) 
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this probability is at most 2[t k — tfc-i] K d 2 ■ Therefore, 

P{for some Xk-i with ||xfc_i|| < 2C\t k -\ and (x k -\,u) > d k -i, 
B (xk-i, tk-i; u, — r) occurs, but £ k) i fails} 



< 2K A [tk-i] d VA[tk-t k -i]- K - d - 2 < K^ h ' 



It turns out that in the estimates for £ kj 2,£k,3 and £k,5 we can ignore the type 
of the particle at (xk-i, tk-i)', we just need that this space-time point is occupied. 
For £k,2 we again shift by (—Xk-i,—tk), sum over the possible values of Xk-i and 
apply (3.7). This gives 

P{for some Xk-i with ||xfc_i|| < 2C\t k -\, (xk-i,tk-i) is occupied but £ kj 2 fails} 
< tf 4 [*fc-i] 0-) = for all y E (C 2 /4)(£ fc - t k -i)u + C([logt k ] 2 )} 

for large r, because the Na(u, 0— ) are independent. 

Next, for £ k ^ we use that on T> k _i, the distance between x k -\-\-C{2C\{t k — t k -i)) 
and the complement of S(u, — r) is at least 

(x k -i,u) + r - 2\fdC\(t k - tk-i) > d k -i + r- 2\ / dC 1 (t k - t fc -i) 

1 , 

Thus, if we take the restriction (xk-i,u) > dk-i into account we find that 

P{for some Xk-i with ||xfc_i|| < 2C\t k -\ and (xk-i,u) > dk-i, 
x k -i is occupied at time t k -i but £ k) s fails} 

< £ £ E P{S tk ^=z-y} 

<.2C\t k 
(x,u)>d k . 



\\x\\<2C 1 t k ^ 1 y^S(u-r) z6x+c ( 2Ci(tfe _ tfc _ i) ) 



< Kstt.lt, - t^PqSt^ || > id fc _! + r} 

< Kftjfexp \ - K 8 + (dk ~ 1 ^ r) " ] (by (2.42) in [KSa]) 

< f** -2 . (3-24) 

The estimate for £^4 comes from Theorem A, or rather Theorem 1 in [KSb], 
which is the basis for the right hand inclusion in Theorem A. Indeed, we have, 
again by summing over the possible values of x k -\ and using (3.9), 

P{for some x k -\ with ||xfc_i|| < 2C\t k _\, (x k -i,t k _i) is occupied but £ k ^ fails} 
< K±\t k -i\ /^P or {there are B-particles outside C(Ci(t k — £fc-i)) 

at some time <t k — t k -i}. 

(3.25) 
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To estimate the probability on the right we argue as in the proof of Theorem 3 
of [KSb]. If a particle has type B at some time s < tk — tk-i and is outside the 
cube C{C\{tk — tk-i)) at that time, then by symmetry of the random walk {S}, 
the particle has a conditional probability given T s at least 1/2 of being outside 
C(Ci(tk — tk-i)) at time tk — tk-i- Therefore (with E or denoting expectation with 
respect to P or ), 

E or {number of particles outside C{C\(tk — tk-ij) at some time 

during [0, t k - t k -i}} 
< 2£ ,or {number of particles outside C(Ci(tk — tk-i)) at time tk — tk-i}- 

The expectation in the right hand side here is exponentially small in (tk — tk-i) by 
Theorem 1 of [KSb] and is an upper bound for the probability in the right hand of 
(3.25). Thus the left hand side of (3.25) is at most % K ~ 2 again. 

The probability involving 5 is also 0(t^ K ~ 2 ). This can be shown by large 
deviation estimates for random walks, analogously to the terms involving £jr 3 . 

This provides the necessary bounds of the summands in (3.23). Finally, we have 

P{V Q fails} < P{N A (x,0-) = for all x with < if 3 logr} 

= exp [ - /j A K 9 [K 3 log r] d ]. (3.26) 

Thus we can take K 3 = K^(K, d, /ja) so large that the left hand side of (3.26) is at 
most r - K ~ 1 for all r > 3. (3.20), (3.26) and (3.23) together now show that the left 
hand side of (3.16) is for large r at most 

oo oo 

k=i k=i 

For any vector v G M d , we define 

v ± = v ± (u) := v — (v,u)u. 

We further introduce the following (semi-infinite) cylinders with axis in the direction 
of it, for a, (3 e R and 7 G M, d a vector orthogonal to u (see Figure 1): 

r(a, /3, 7) = r(a, /3, 7, u) := {x G R d : (x, u) > a, HaH" - 7|| < 0) 

and the events 

Q(a, (3, 7, V, t) = Q(a, (3, 7, V, t, u) 

: = {in the process V there is a S-particle in T(a, (3, 7) at time t}. 
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Figure 1. The shaded region represents a cylinder T(a, (3, 7); it extends 
to infinity on the upper right. 

The last definition will be used with V taken equal to some half-space, full-space 
or the original process. 

We remind the reader that V h (u, — r), the (w, — r) half-space process, only uses 
particles which at time are in the half-space S(u, — r) = {x : (x,u) > — r}. We 
shall work a great deal with the process T )h {u : —C§k(s)), where 

K{s) = [{s + l)\0g{s + l)] 1 / 2 

and C 5 is a constant to be determined below (see the line preceding (3.63)). We 
make several more definitions: 

h*(s, u) : = h(s, u, — Csk(s)) = max{(i, u) : B h (x, s; u, — Csk(s)) occurs} 
= max{(i,«) : x is occupied by a S-particle in 

V h (u, -C 5 k(s)) at time s}. 

(3.27) 
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We order the sites of Z d in some deterministic way, say lexicographically, and take 

£*(s, u) := the first x in this order for which B h (x, s; u, — C^,k(s)) occurs 

and which satisfies (x,u) = h*(s,u). (3.28) 

Thus, h*(s,u) is the furthest displacement in the direction of u among the B- 
particles in the process V h (u, — Csk(s)) at time s and £*(s, u) is the first site occu- 
pied by a .B-particle in this process at time s on which this maximal displacement 
is reached. We shall write m*(s,u) for [Z*(s, it)]- 1 so that we have the orthogonal 
decomposition 

t(s,u) = h*(s,u)u + m*(s,u). (3.29) 

The following proposition contains our principal "subadditivity" property. If we 
take (3 = oo, that is, if we only look at its statement about displacements in the 
direction of u, then the proposition says that (up to error terms) the maximal dis- 
placement in the direction u at time s + 1 + C§K(t) in the process V h {u, —C^k^s + 
t + CQK,(t))) is stochastically larger than the sum of two independent such displace- 
ments, which are distributed like the maximal displacement in V (u, — Csk,(s)) at 
time s and the maximal displacement in V h {u,—C^K(t)) at time t, respectively 
(see Corollary 5 for more details). The basic idea of the proof is that if £* is a 
point where V h (u, — Csk(s)) achieves its maximum displacement in the direction 
u at time s, then we can start a new half-space process "near" £* which is nearly 
a copy of V (u, — C§n{t)\ and which is nearly independent of the first process 
T >h {u : —C$k(s)). If we run the second process for t units of time the sum of the 
displacements in the direction of u in the first and second process is nearly a lower 
bound for the displacement of the original process at time s + t. 

Proposition 3. Let u E 5' d ~ 1 ,a, / 9 > and 7 s ,7t G R d orthogonal to u. For any 
K > there exist constants < C5 — C%, sq < 00, which depend on K , but are 
independent of u £ 5' d_1 and of a, ft, 7 S , 7 t , such that for 

s < s < t and tlogt < C 7 s 2 (3.30) 

it holds 

P{Q (a, f3, 7s + 7t , V h (u, -C 5 k(s + t + C 6K (t))) ,s + t + C 6K (t)) } 
> / P{h*(s, u) E dh, m*(s, u) — 7 S G dm} 

J0<h<oo JmeR d 

P{Q(a -h,/3-d,>yt-m, V h (u, -C 5 k(*)) , t) } 



'0<h<oo Jm 

X 

-K-l 



-C 8 s~ 

(3.31) 



Proof. The constants C\ and so will be fixed later on. K{ will be used to denote 
other auxiliary constants. For the time being we only do manipulations which do 
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not depend on the specific value of the Ci,Ki. We break the proof up into three 
steps. The first two steps reduce the proof of (3.31) to estimating the probabilities 
of a number of small events. These estimates will be carried out in step 3. 
Step 1. Run V h (u, -C 5 k(s)) till time s. Let h*(s,u) = h G R, t (s,u) = y G Z d . 
Set y := \_y+AC*,K{t)u\ (the meaning of this last notation is that we take the largest 
integer for each coordinate separately). Next we run the (u, (y,u) + 2C 5 k(£)) half- 
space process starting at the space-time point (y, s + Cqk{£)^ for t units of time. 
This latter half-space process will be shown to be almost an independent copy of 
the translate by (y, s + Cqk{€)) of V h {u, —C$k{$)). Define z(s, t) to be the nearest 
site in Z d to y which is occupied at time s + C 6 K(t) by a particle which started at 
time in the halfspace S(u, (y, u) + 2C 5 k(£)). It will be useful to define for general 
v G Z d 

z v = the nearest site on Z d to v := [v + 4C^,K(t)u\ which is 
occupied at time s + Cqk{€) by a particle which started 
at time in S(u, (v,u) + 2C 5 «(t)). (3.32) 

Thus, z v has the same relation to v as z(s,t) has to y. In particular, z y = z(s,t). 
We can now define, still for any v G Z d , the sets 

A\(s, t) = {x : x is occupied by one or more S-particles at time 
s + t + CgK(t) in the process 

A2(s,t, v ) = {x : x is occupied by one or more S-particles at time 

s + t + Cen(t) in the (it, (v, u) + 2C^n{t)) half-space 
process starting at (v, s + C 6 K(t)}, (3.33) 

and 

A s (t) = {x : x is occupied by one or more S-particles at time t in 
an independent copy of the process V h (u, —C 5 n(t))}. 

We stress that A 3 is defined by means of a new copy of all initial data and paths. 
It is independent of the processes we worked with so far. 

Our aim is to prove the following two statements, and to show that they imply 
(3.31). The first statement is that outside an event of probability at most s~ K ~ 1 
it is the case that 

Ai{s,t) DA 2 (s,t,y). (3.34) 

The second statment is that 

Ai(s, t) — y is at least as large as A 3 (t), outside 

an event of probability at most s~ K ~ 1 (3.35) 
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(still y = £*(s, u) in these relations). The relation (3.35) is stated somewhat impre- 
cisely, but a precise version will be given below (see (3.51)). In this first step we 
shall reduce the proofs of (3.34) and (3.35) to a number of probability estimates. 

To begin with the inclusion (3.34), we claim that this holds on the intersection 
of the event 

{{y, u) > 0} fl {z(s,t) is occupied by a S-particle at time 

s + C 6 K(t) in V h (u, -C 5 k(s))} (3.36) 

with the event (see (3.4) for x ) 

{\\x \\<K 3 logs}. (3.37) 

This follows from two applications of the monotonicity property in Lemma C. In- 
deed, under (3.37) (and s > s± for a large enough si), both the (-u, — C5«:(s)) and 
the (w, — C$k(s + t + CsK(t))) half-space process begin with B-particles at xq. One 
application of the monotonicity property therefore gives us that (under (3.37)) 
V h (u, -C 5 K(t + s + C & K{t))) has more S-particles than V h (w, — C$k>(s}) at each 
space-time point, and therefore 

Ai(s,t) D {x : x is occupied by one or more -B-particles at time 
s + t + C 6 K(t) in the process V h (u, — C 5 k(s)) }. 

(3.38) 

For the second application we recall that z(s, t) is occupied at time s + C§K(i) by a 
particle which started in S(u, (y, u) + 2Csk(£)), and in fact is the closest occupied 
site to y with this property. To run the (u, (y,u) + 2C§K{t)) half-space process 
starting at (y, s + CQK,(t)\ and to find ^(s, t, y) we first remove all particles which 
at time were in the halfspace {x : (x,u) < (y,u) + 2C^n{t)}. After that we reset 
to A the types of all particles not at z(s, t) at time s + Cqk{£) and give all particles 
at z(s, t) type B. Note that in the first step all particles which do not belong to 
V h (u, —C 5 k(s)) are removed, since 

-C 5 K(t) < 2C 5 K(t) < (y,u) + 2C 5 K(t) (on (3.36)). 

Thus, (at time s + C^K^t)) after both steps, all remaining particles are also in 
V h {u, — C5«(s)), and the particles which have type 5, i.e., only the particles at 
z(s,t), also have type B in V h (u, — C$k(s)) (still on the event (3.36)). By virtue 
of the monotonicity property of Lemma C, at time s + t + CQK,(t), any S-particle 
present in the (u, (y, it) + 2C 5 n(t)) half-space process starting at (y, s + C 6 K(t)) is 
also a B-particle in V h (u, —C 5 k(s)). Therefore, on the event (3.36), 

A2(s,t, y) C {x : x is occupied by one or more -B-particles at 

time s + t + CeK,(t) in the process V {u, — Csk(s))}. 

(3.39) 
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Combining (3.38) and (3.39) gives (3.34) on the intersection of the events (3.36) 
and (3.37). We postpone the proof that this intersection indeed has probability at 
least 1 — s~ K ~ 1 to step 3. 

To prepare for the desired precise form of (3.35) we shall prove that there exist 
constants K\ and s 2 such that for t > s > s 2 , A any non-random subset of Z d , and 
any fixed v G Z d , 



P{A 2 (s,t,v) intersects A} > P{(v + A 3 (t)) intersects A} - K 1 t~ K ~ d ~ 1 . (3.40) 



To prove this inequality we remind the reader that A 2 (s,t, v) is the collection of 
sites where S-particles are present at time s + t + CeK(t), if one starts at time 
s + C§K(t) in the state obtained by removing the particles which started outside 
S(u, (v,u) + 2CsK,(t)) at time 0, and by resetting all particles not at z v to type 
A, while setting the type of the particles at z v to B. To find the distribution of 
A 2 (s, t, v) we must first describe the state at time s + C§K(t) (after the removal of 
particles and resetting of types) in some more detail. First let us look how many 
particles there are at the various sites, irrespective of their type. We began at 
time with Na(w,0— ) particles at w, for w G S(u, (v,u) + 2C^K{t)) and with 
particles at any w outside S(u, (v,u) + 2C$K(t)). The A^4(m;,0— ) are i.i.d. mean 
Ha Poisson random variables. We let these particles perform their random walks 
till time s + CgK(t). Let us write A^(w, s + C6K(t)) for the number of particles (of 
either type) at w at this time. By properties of the Poisson distribution, all the 
N(y + w,s + C 6 K(t)), w G Z d , are still independent Poisson variables, but 

EN(v + w,s + C 6 K(t)) 

= VAP{S a +c e K(t)=v + w-v/}=:u(y,w,s,t). 

w 'es(u,{v,u)+2C 5 K(t)) ( 3 - 41 ) 



z v is now the nearest lattice point to v which is occupied by some particle at 
time s + Cen(t). We then reset all particles not at z v to type A, and the ones at 
z v to type B. If we shift everything by ( — v, — s — Cqk{€)) (i.e., move (w,r) to 

(w — v, r — s — C6«(t))), then, at (w, 0) we have M(w) := A^(u + w, s + C§K(t)) 
particles, all of which will be reset to type A, except those at too := the nearest 
lattice site to the origin with M(w) > 0. In fact, wq = z v — v. The M(w) are 
independent Poisson variables, and M(w) has mean v(v, w, s, t). It follows from the 
definition of A 2 (s, t, v) and of the (u, (v, u) + 2C 5 n(t)) half-space process started at 
(y, s + C 6 K(t)) that A 2 (s, t, v ) — v has the distribution of 

{x : there is a S-particle at x at time t in this shifted system}. (3.42) 

For the w G S(u, —C^K{t)) : the means u(v, w, s, t) are close to \ia- In fact, it follows 
from (3.41) that for w G S(u, —C^K^t)) and t > to V s, for some to (independent of 
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v,w,u), 

fJ,A>v(v,W,S,t)>fJ, A l- Yl P { S s+C 6 K(t) = 

w:{w,u)>CsK(t) — d 

>Ha[1- P{(S 8+ c 6K(t ),u) > C 5K (t) - d}] 

> n A [l -K 2 exp[-K 3 Cllogt}]. (3.43) 

for some constants K 2 , K 3 that depend on d, D only (see (2.42) in [KSa] and recall 
that we assume (3.30)). From now on we take C 5 so large that for large t 

(jl a [1 -t~ K - 2d - 1 ] < v(v,w,s,t) < /I A for all w G S(u, — C 5 «(*))- (3-44) 

It suffices for this that K 3 C\ > K + 2d + 2. We may have to raise C5 in the proof 
of (3.62) and (3.64) in step 3, but that can only improve the present estimates. 
With such a choice of C 5 the distribution of the particle numbers {M(w) : w G 
S(u, —C$K,(t)) nC(3Cit)} differs in total variation from the distribution of an i.i.d. 
collection of mean \ia Poisson variables on S(u, — C$K{i)) fl C(3Cit) by at most 

fi A t- K - 2d - 1 < K 4 t- K - d ~ 1 (3.45) 

io€s(u,-C 5 K(t))n C(3dt) 

for some constant K4 = K^^a, d). 

Now consider an auxiliary process which starts at time with Na(w 7 0—) par- 
ticles only at the vertices w G S(u, — CsK(t)) fl C(SC\t), and with no particles 
outside this set. Let too be the nearest vertex in S(u, — C5«(t)) to the origin 
with iV^iuo, 0— ) > 0. Take the type of all particles not at wq equal to A and 
the type of the particles at wq equal to B. If u>o lies outside S(u, — C^K^t)) fl 
C(3Cit), then this auxiliary process has never any S-particles. On the other 
hand, if wo G S(u, —CsK(t)) flC(3Cit), then the auxiliary process is obtained from 
V h {u, —C 5 n(t)) by removing at time all particles in S(u, — C 5 n(t)) \C(SCit). 
Finally, let 

A±{t) = {x : there is a S-particle at x at time t in this auxiliary system}. 

From our considerations above (in particular (3.42), (3.45)) we have that 

P{A 2 (s,t,v) intersects A} > P{v + A 4 (t) intersects A} - K 4 t~ K ~ d ~ 1 . (3.46) 

Indeed, were it not for the fact that Na(w,0— ) is a Poisson variable of mean /j,a 
instead of v{v, w, s,t), the auxiliary system would be obtained from the system in 
which A 2 (s, t, v) is computed by translation by ( — v, — s — 6*6 and by removing 
the particles outside <S(tt, — C$K{i)) flC(3Cit). The term — K 4 t~ K ~ d ~ 1 corrects for 
increasing the mean from v(v, w, s, t) to \ia- 



SHAPE THEOREM FOR SPREAD OF AN INFECTION 



23 



To come to (3.40) we still want to prove the inequality 

P{v + A 4 (t) intersects A} > P{v + A 3 (t) intersects A} - K$t~ K ~ d ~ l . 



(3.47) 



This follows from the fact that if, in V h (u,—C§K{i)) all S-particles stay inside 
C{2C\t) during [0,t], and no particle which starts outside C{3C\t) at time enters 
C{2C\t) during [0, t], then the particles which start outside C{3C\t) do not interact 
with any particle, and do not cause the creation of any S-particles during [0, t] 
(compare the argument for (2.36) in [KSb]). In these circumstances V h (u, —CsK{tj) 
has no more S-particles at time t than the auxiliary process, which is obtained by 
removing the particles which start outside C(3C\t) at time 0, as described above. 
Therefore 

\P{v + A 4 (t) intersects A} - P{v + A 3 (t) intersects A| 
< P{w i S{u, -C 5 K(t)) n C(3Cit)} 

+ P{some S-particles in V h (u, -C 5 «(t)) leave C{2C\t) during [0,t]} 

some particles which start outside C{3C\t) 

enter C{2C x t) during [0,t]}. 

(3.48) 

The first term in the right hand side here is trivially o(t~ K ~ d ~ 1 ) (compare (3.26)). 

To estimate the second term in the right hand side of (3.48) we first observe 
that if xo, the nearest occupied point to the origin at time 0, lies in the half-space 
S(u, — C 5 n(t)), then at all times the S-particles in V h (u, —C 5 n(t)) are also B- 
particles in (see (3.3)). Therefore a decomposition with respect to the value of 
xo and an application of (3.7) show that the second term in the right hand side of 
(3.48) is (for t > some tz) bounded by 

-P{||^o|| > K± logt} + ^ ^aP° t {there are S-particles 

||x||<K 4 logt 

outside C(dt) during [0,t]} 
< K G r K - d - 1 (3.49) 

(see (3.3), (3.26) and the estimates for (3.25)). 

The third term in the right hand side of (3.48) is at most 

E{N A (w, 0-)}P{sup ||5 r || > ||w|| - 2dt} 

w$C(3C\t) r - f 

< 8dfiAex.p[-K 7 \\w\\] < Kst^^xpl-Kgdt] 

w£C(3Cit) (3.50) 



(see (2.42) in [KSa]). Thus (3.47) and (3.40) hold. 
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Step 2. We wish to prove the following precise version of (3.35): for t > s > 
so,tlogt < C7S 2 and for some constant -fTio, independent of s,t,u, 

P{Ai(s,t) intersects A} 

> P{ (£*(s,u) + A 3 (t)) intersects A} - K w t~ K ~ 1 . 

(3.51) 

To this end we define the following events for any vector in v E Z d : 

T\{v) := j during [0, s] in the process V h (u, — C 5 k(s)) all the S-particles 
stay in the set C(2Cis) fl {x : (x, u) < (v, u) + C 5 k(£)} j, 

1 2 {v) := I none of the particles which were at time in the half-space 
S(u, (v, u) + 2C 5 n(t)) = {x : (x, u) > (v, u) + 2C 5 K(t)} enters 
the set C{2C 1 s) n {x : (x,u) < (v,u) + C 5 K,(t)} during [0, s]}. 

The following independence property is crucial for our argument: Let J(v) be an 
event which depends only on v E Z d and the particles which start in <S(tt, (v, u) + 
2C 5 n(t)) at time 0, and the paths of these particles. Then 

P{t(s,u) = v,Xi(v) t l2(v) t J(v)} 

= P{e*{s,u) = v,l 1 {v),l 2 {v)}P{J{v)\l 2 {v)}. (3.52) 

The important feature here is that in the last conditional probability v is a constant, 
without relation to £*(s, u). To see (3.52) we note that on the event X\ fl X 2 none 
of the particles which start in S(u, (v,u) + 2C^K,(t)^ coincides with any S-particle 
during [0, s] . Therefore, changing the paths of any of the particles which start in 
S(u, (v,u) + 2C$K(t)) has no influence on the types of any of the other particles 
during [0, s] (and of course no influence on the paths of these other particles), as 
long as we stay on X\ flX 2 (compare the argument for (2.36) in [KSb]). In particular, 

P{t(s,u)=v,X 1 (v)\X 2 (v),J(v)} = P{t(s,u)=v,X 1 (v)\X 2 (v)}. 

This is clearly equivalent to (3.52). 
We take 

J{v) = {A 2 (s, t,v) intersects A}, 

where A is some non-random set in Z d . By definition, A 2 (s,t,v) depends only 
on v and the particles which start in the half-space S(u, (v, u) + 2C 5 K(t)) . Thus 
also Jiy) depends only on v and this last collection of particles and their paths. 
(This is true despite the fact that we talk about S-particles in the definition (3.33). 
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Indeed, these are B-particles in (ix, (v,u) + 2C^n{t)) half-space process, started at 
(y, s + Cqk,{1:)\ and the types of these particles are reset at time s + CQK(t) and 
after that do not depend on particles which started outside <S(ix, (v, u) + 2C5 «(£)).) 
With this choice of J we obtain from (3.52) for every fixed v 

P{£*(s,u) = v,Zi(v),I 2 (v), A 2 (s,t, v) intersects A} 

> P{£*(s,u) = v,l 1 (v),l 2 (v)}[P{A 2 (s,t,v) intersects A} - P{l 2 (v)}] + . 

(3.53) 

We shall show in step 3 that for suitable choice of constants < Ki = Ki(K, d) < 
oo, independent of s,u and v, it is the case that for the process V (u, —C5k(s)) 

P{t(s, u) e C(2C lS ), 2J(r (s, u))} < K llS - K -\ (3.54) 

P{l c 2 (v)} < K llS - K - d -\ (3.55) 

and 

P{(3.36) (with y = t(s,u)) fails or (3.37) fails} < K 11 s~ K - 1 . (3.56) 

In the remainder of this step we only show how to complete the proof of (3.51) and 
the proposition from the estimates (3.54)-(3.56). To this end we apply (3.53). By 
using (3.53), (3.55) and (3.40) we get 

P{t(s,u) = v,A 2 (s,t,t(s,u)) intersects A} 

> P{£*(s,u) = v,Ii(v),I 2 (v), A 2 (s,t,v) intersects A} 

> P{t(s,u) =v} - P{t(s,u) =v,lf(v)} - P{t(s,u) =v,ll(v)} 

x [P{A 2 (s,t,v) intersects A} - P{l%(v)}] + 

> P{£*(s,u) = v}P{A 2 (s,t,v) intersects A} 

- P{t(s,u) = u,2j(u)} - 2K 11 s~ K - d - 1 

> P{£*(s,u) = v}P{v + A 3 (t) intersects A} 

- P{t(s,u) = v,lZ(v)} - (2K U +K 4 + K 5 )s- K - d - 1 . 

(3.57) 

Now recall that Ai(s,t) D A 2 (s,t,£*(s,u)) on the intersection of (3.36) and (3.37). 
Summing (3.57) over all v £ C(2Cis), and using (3.54) and (3.56), therefore gives 

P{Ai(s,t) intersects A and (3.36), (3.37) occur} 

> P{A 2 (s,t,t(s,u)) intersects A} - P{(3.36) or (3.37) fails} 

> P{Z*(s,u) = v}P{v + A 3 (t) intersects A} -K 12 s~ K - 1 . 
veC(2C lS ) (3.58) 
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Finally, since £*(s,u) is the location of a S-particle at time s in V h (u, —Csk(s)) : 
we have, essentially as in the estimate for P{£^ 4 } in (3.25) and the lines following 
it, or the estimate of the second term in the right hand side of (3.48) 

P{t(s, u) i C(2C lS )} < P{\\x \\ > ds A C 5 k(s)/VcI} 

+ [i-AP or {there is a 5-particle outside C{C\s) at time s} < Ki 3 s~ K ~ 1 . 

\\x\\<c lS (3.59) 

Consequently 

P{Ai(s,t) intersects A} > ^ p {^i s -> u ) = v}P{v + A 3 (t) intersects A} 

vei d 

-(K 12 + K 13 )s- K -\ 

(3.60) 

This is the desired (3.51). 

(3.31) is just the special case of (3.60) with A = F(a,P, r y s + 7 t ). Indeed, 
{Ai(s,t) intersects A} is the event that there is a S-particle in A at time s + 
t + C 6 K(t) in the process V h (u, -C 5 k(s + t + C 6 «(t))) . For A = T(a, /3, 7 S + 7t) this 
event is also denoted by Q(a, (3, 7 s + 7t, V h (u, — C 5 K(s + t + C 6 «:(t))), s + t + C 6 K(t)). 
Thus, the left hand sides of (3.31) and (3.60) are the same for this choice of A. We 
leave it to the reader to check that the right hand side of (3.60) is at least as large 
as the right hand side of (3.31), provided we choose C$ > Kyi + K13. 

Step 3. Here we prove the relations (3.54)-(3.56). Note that (3.56) also supplies 
the missing estimates for (3.34), to wit, P{(3.36) and (3.37) hold} > 1 — s~ K ~ 1 . 
Now we start on (3.54). First 

P{in V h (u, -C 5 k(s)) some B-particle leaves C{2ds) during[0, s]} = 0{s~ K - d - 1 ) 

(3.61) 

(see (3.49)). In addition, by definition of l*(s,u), (l*(s,u),u) = h(s, u, — Csk(s)) . 
Thus 

P{in V h (u, —C 5 k(s)), during [0, s] all S-particles stay in C(2C\s), but some 

of them leave {x : (x, u) < (£*(s, u),u) + Csk(s)} } 
< P{in V h (u, — C5«(s)), at some time r < s there are 
B-particles at some v G C(2Cis) with 

(v,u) > max.{{x,u) : there is a S-particle at x at time s} + Csk(s)}. 

(3.62) 

This last event can happen only if some S-particle reaches a vertex v G C(2C\s) 
before time s and then this particle moves to some x at time s with (x,u) < (v, u) — 
C 5 k(s). The probability that such particle started outside C(3Cis) is bounded by 
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the third term in the right hand side of (3.48), with t replaced by s. Therefore, the 
right hand side of (3.62) is at most 

(third term in right hand side of (3.48) with t replaced by s) 

+ Yl E{N A (w,0-)P{ sup \\S ri - S r2 \\>C 5 K(s)/Vd} 
w eC(3C lS ) o<n,r 2 < s 

< Itss*- 1 exp[-K 9 C lS ] + K 14 (3C 1 s) d exp[-K 15 C%logs], (3.63) 



by (3.50) and by (2.42) in [KSa]. Together with (3.61) this proves that we can take 
C5 so large that (3.54) holds. As observed after (3.44) we can even choose C5 so 
that (3.44) is also valid. Once we have chosen C5 we fix 



C 6 = (3.64) 



As for (3.55), we have 



P{n(v)}< E E{N A (w,0-} 

weS(u,(v,u)+2C 5 K(s)) 

x P{sup ||5 r || > C 5 K(t) V (|H| - 2C 1 s)}. 

r ^ s (3.65) 

We leave it to the reader to show that this is 0(s~ K ~ d ~ 1 ) for t > s and large 
enough C 5 (again by (2.42) in [KSa]). 

Finally, to prove (3.56), we note first that P{(3.37) fails} = 0(s~ K ~ 1 ), provided 
K 3 = K 3 ((xa, d) is taken large enough, just as in (3.26). Next, 

P{(£*(s,u),u) < 0} < P{h(s,u,-C 5K (s)) < C 4 s} 

< [C 5K (s)]- 2K - 2 < s~ K - 1 (3.66) 

for large s, by virtue of Lemma 2 with K replaced by 2K + 2. Lastly, we have to 
show that for the choice of C% in (3.64) 

P{z(s, t) is not occupied by a S-particle in V h (u, —C 5 k(s)) at time s + C 6 n(t)} 
< P{z(s,t) is not occupied by a S-particle in V h (u, —C 5 k(s)) at time 

s + C 6 K(t), but z(s,t) E t(s,u)+C{C 2 C 6 K(t)/2} 
+ P{z(s,t) <£e*(s,u)+C(C 2 CeK(t)/2} 
= 0(s- K - 1 ). (3.67) 

The first inequality here is obvious. The bound 0(s~ K ~ 1 ) for the middle member of 
(3.67) is formulated as a separate lemma, because the same argument will be needed 
once more in the next section. To see that (3.67) follows from Lemma 4 below, recall 
that z(s, t) is occupied at time s + C§n{t) by some particle which started at time 
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in S(u, (£*(s,u),u) + 2C$K(i)) (by definition of z(s,t)). In particular there is 
some particle at z(s, t) at time s + CQK.it) , so that z(s, t) is occupied in V? at time 
s + Cqk{€). Also, £*(s, it) is occupied by at least one B-particle in V (u, — Csk(s)) 
at time s. So Lemma 4 with s = s + C§K(t) and y{s) = £*(s,u) (and Cq as in 
(3.64)) shows that the middle member of (3.67) is at most 

P{£*(s,u) i C(2C lS )} + P{(£*(s,u),u)) < C 4 s/2} + s~ K_1 
+ P{z(s,t) t t(s,u)+C(C 2 C 6 K(t)/2)}. 

(3.68) 

Note that we used the second part of condition (3.30) here; we have to choose CV 
small enough to make sure that (3.72) holds for s — s = C 6 K(t). The first two terms 
in (3.68) are 0(s~ K - 1 ), by virtue of (3.61) and (3.66). The fourth term is bounded 
by 

P{z(s,t) i e*(s,u)+C(C 2 C 6 K(t)/2)} < P{\\z(s,t)-£*(s,u)\\ > AC b n{t) - 1} 
(because C 2 C 6 /2 = 8C 5 and \\£*(s, u) - £*(s, u) \\ < 4C 5 n(t)) + 1) 
< P{£* (s, u) i C(2C lS )} + P{t(s, u) e C(2C lS ), and none of the sites 
in £*(s, u) + C{4C§K(t) — l) are occupied at time s + C^nit) by 
a particle which started in S{u, (£*(s,u),u) + 2C 5 K,(t))}. (3.69) 

We already saw in (3.59) that the first term in the right hand side is 0(s~ K ~ 1 ). As 
for the second term in the right hand side, this is by a decomposition with respect 
to the possible values of £*(s,u), analogously to (3.9), at most 

if i 6 ^ P{none of the sites in v + C{AC^,K{t) — l) is occupied at time 

veC(2C lS ) 

s + CQn(t) by a particle which started in <S(-u, (v, u) + 2C^K(t))}. 

(3.70) 

However, the numbers of particles at sites v + w at time s + Ce,K(t) which started 
in S(u, (v, u) + 2C$K{i)) are independent Poisson variables with means u(v, w, s, t) 
given in (3.41). By the estimate (3.43) we have u(v,w, s,t) > for (w,u) > 

and all v (and t large enough). Therefore (3.70) is at most Kns d exp[— Ki%K d {t)iJiA\- 
This proves the bound 0{s~ K ~ 1 ) in (3.67), and therefore (3.56) is reduced to the 
next lemma. 

Roughly speaking, the next lemma guarantees that if a certain vertex y(s) has a 
5-particle in the half-space process V h (u 1 —C^nis)) at a time s, then a little later 
all occupied sites "near" y(s) will actually have a S-particle in V (u, — C5«;(s)). 

Lemma 4. Let s,s and t be such that 

s < t (3.71) 
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and 

^n{t)<s-s<^-s. (3.72) 

Let u E >S d_1 be fixed and let y(s) E Z d be a random point (that is, y(s) may depend 
on the sample point a). Define the event JC(y) by 

JC(y) := {there exists a site z E y + C(C2(s — s)/2j such that 
at time s, z is occupied in , but is not occupied 
by a B-particle in V h (u, — C$k(s))}. (3.73) 

Then for each K > there exists an si = si(K) (independent of u) such that 

P{B h (y(s),s;u, -C 5 K(s))nlC(y(s))} 

< P{y(s) i C(2C lS )} + P{(y(s),u) < ±C±s} + s"^ 1 (3.74) 

fors> si (see (3.13) for B h ) . 

Proof. Assume that the space-time point (y, s) is occupied by some particle in 
V h (u, — C5«(s)). We can then define the following auxiliary events: 

lCi(y) :— {there exists a site z E y + C(C2(s — s)/2) such that (z,s) 

is occupied in V, but is not occupied in V (u, — C5«(s))}, 
^2(2/) '■— {there exists a site z E y + C(C 2 (s — s)/2) such that (z,s) 

is occupied by an A-particle in the full-space process 

starting at (y, s)}, 
ICs(y) := {there exists a site z E y + C(C 2 (s — s)/2) such that (z,s) 

is occupied by a -B-particle in the full-space process starting 

at (y, s), but occupied by an A-particle in the (u, —Cr,K(s)) 

half-space process starting at (y, s)}, 
/C 4 (y) := {in the full-space process starting at (y, s) some S-particles 

leave y + C(2Ci(s — s)) during [s,s]}, 
^5(2/) := {some particles which start outside S(u, —Csk(s)) at time 0, 

enter y + C(2Ci(s — s)) during [s,s]}. 

We shall first show that 

3 

B h {y, s; u, -C 5 «(s)) n K{y) c |J £*(?/) and /C 3 (y) C /C 4 (j/) U £ 5 (y), (3.75) 

i=l 

and then estimate P{y(s) G C(2Cis), (y(s), u) > C 4 s/2, IC t (y(s))} for 1 < i < 5. To 
prove the first part of (3.75) , consider a sample point for which B h [y,s;u, — C§k{s)) D 
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JC(y) occurs and let z be a site in y + C(C2(s — s)/2) such that (z, s) is occupied in 
V* , but is not occupied by a S-particle in V h {u, —C$k{s)). Then it may be that 
(z,s) is not occupied at all in V (u, —C^^s)). This would mean that /Ci (y) occurs. 
If this fails, then (z,s) is occupied in V (u, — 65^(5)), necessarily by an A-particle. 
We claim that (-2,5) is then also occupied by an A-particle in the (w, — Csk(s)) 
half-space process starting at (y, s). This is so, because starting at (y, s) does not 
remove any particles, but it may change some types. But on B h (y, s;u, —C^k^s)), 
y has already at least one S-particle at time s in V h (u, —C^^s)). Thus the re- 
setting at time s only changes some types from B to A, and since z already has 
type A at time Is in V h {u, —C§k(s)), it will (by Lemma C) also have type A at 
time s in the (u, — C^k{s)) half-space process started at (y, s), as claimed, (z, 3) is 
also occupied in the full-space process starting at (y, s) (since it is occupied in the 
full-space process, starting at (0,0)). The type at (z, s) in this process may be A, 
in which case £2(2/) occurs, or B, in which case ICs(y) occurs. This proves the first 
inclusion in (3.75). 

The second part of (3.75) follows from the argument given for (3.47). ICs(y) 
requires that there are particles in y + C(C 2 (s — s)/2) which have different types 
at time sin the full-space and in the (u, — Csk(s)) half-space process, both starting 
at {y,s). This means that in the full-space process starting at (y, s) the type of 
some particle which is in y + C(C2(s — s)/2j at time 's is influenced by particles 
which started outside S(u, — Csk(s)) at time 0. However, this can happen only if 
in the full-space process starting at (y, s), these particles meet some B-particles 
during [s,s\. In turn, this can happen only if )C 4 (y) or jC$(y) occurs. This proves 
the second inclusion in (3.75). 

Our next task is to find bounds for 

P{y(s) G C(2C l8 ), (y(s), u) > C 4 s/2, B h (y(s), s; u, -C 5 k(s)), Ki(y(s))}, i = 1,2, 4, 5. 
For i = 1 we have 

P{y(s) e C(2C lS ), (y(s), u) > C 4 s/2, /d (?/(*))} 

< Yl EN A (w,0-)P{w + S~ s = z} 

w$s(u,-C S K(sj) (z,u)>C f s/2-C 2 (s-s)/2 
z6C((2Ci+C 4 )s) 

< Yl »aP{\\S~ s \\ >C 4 s/(4Vd)}+ J2 »aP{\\Ss\\ > HI/2} 

wec((4C 1 +2C 4 )s) ■wfc((4C 1 +2C 4 )s) 

< s~ K ~ l (3.76) 

for all s > some s\ = si(K). In the first inequality we used that ||z|| < \\z — y(s)\\ + 
\\y(s)\\ <C 2 (s- s)/2 + 2ds < (2Ci + C 4 )s, by virtue of (3.72) and the inequality 
C2 < C\. In the second inequality we used that for the summands here we have 
1 1 to — z\\ > ((z — w),u)/Vd> [C 4 s/2 - C 2 (s- s)/2 + C 5 K(t)]/Vd > C 4 s/(A^d). 
For the third inequality we use s < (l + C 4 /(8Ci))s plus (2.42) in [KSa]; compare 
(3.24). 
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Next, we remind the reader that P or is the probability measure governing the 
original model, in which one B-particle is added at the origin at time 0. In this 
notation we have, by (3.9) and (2.4), 

P{y(s) e C(2C l s),B h {y(s),s- u, -C7 5 «(a)), /C 2 (y(a))} 

< K 19 s d P or {there exists a ze C(C 2 (s- s)/2) which is 

occupied by an A-particle at time s — s} 

< 2s~ K - 1 . (3.77) 

Again by (3.9) 

P{y(s)eC(2C lS ),)C4y(s))} 

< Ki 9 s d P or {some S-particles leave C{2C 1 (s - s)) during [0,s - s]} 

< s~ K - 1 (by (3.49)). (3.78) 

Finally, 

P{y(s) e C(2C lS ), (y(s), u) > C 4 s/2, IC 5 (y(s))} 

< Yl en a(w, 0-)P{w + S r = v for some r <s} 

w.(w,u)<-C 5 k(s) veC ^ 2Ci+Ci)s ^ 
{v,u)>C 4 s/4: 

<fi A p { SU P \\Sr\\>C 4 s/(4Vd)} 

we c((4C 1+ 2C 4 )s) r<(l+C 4 /(8 Cl ))s 

+ HA Yl P { SU P ll^rll > |k/2||} 

W e((4C 1+ 2C 4 ) S ) -<(l+^4/(8C 1 )) S 

< s~ K - 1 (compare (3.76)). (3.79) 

Together with (3.75) these estimates prove (3.74). ■ 

Corollary 5. For every unit vector u there exists a constant \{u) G [C 4 ,2\fdCi] 
such that 

lim -h*(t, it) = X(u) almost surely and in LP for all p > 0. (3.80) 

t^oo t 

(t runs through the reals here). Moreover, for eachr\ > there exist an exponentially 
increasing sequence {ri\ < n 2 < ■ ■ ■ } = {n\{ri) < n 2 (r]) < . . .} (independent of u) 
such that 

1<^<1 + V, 3>h 



and such that every e > 0, 



h*(n k ,u) - X(u) > e} < oo. (3.81) 

fc=o Uk 
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Finally, for given rj > and rj' < rj and {rij(r})} corresponding to n, we can choose 
the nj(n') such that for some jo < oo, {nj(rf)} ^ { n j(v) '■ 3 — Jo}- 

Proof. The basis for this proof is (3.31) with (3 = oo. Since T(a, oo, 7, u) = {x E 
R d : (x, u) > a} = S(u 7 a), we have 

G(a, 00, 7, V, t) = {in V, at time t, there is a S-particle at some x with (x, u) > a}. 
In particular 

G(a, 00, 7, -C 5 k(s + t + C 6 K(t)), s + t + C 6 K(t)) 

= {h(s + t + C 6 K(t),u, -C 5 k(s + t + C 6 K(t))) > a} 
= {h*(s + t + C 6 K(t),u) > a}. 

Similarly, 

g {a, 00, 7, V h (u, -C 5 K(t)), t) = {h* (t, u)>a}. 
Thus, (3.31) with (3 = 00 says that, under (3.30), 

P{h*(s + t + C 6 K(t),u) >a}> P{hl(s,u) +h* 2 (t,u) > a} - C 8 s~ K ~\ (3.82) 

where h\{s,u) and h\(t,u), are independent copies of h*(s,u) and h*(t,u), respec- 
tively. 

The Corollary will be derived from this relation by more or less standard subad- 
ditivity techniques. To apply these techniques we first derive some simple properties 
of h*(s,u). The first is the following tail estimate: 

P{h*{s,u) > a} + P{\\m*(s,u)\\ > a} < expf-i^c*] for a > 2y/dds, s > s 3 , 

(3.83) 

where S3 = ss(K) is some constant independent of a,u. The second and third 
property are semi-continuity properties in s, namely 



P{inf h* (s + r, u) - h* (s, u) < -a} 

r<t 



< K 3 s~ K + P{\\ sup5 r || >a}< K 3 s~ K + 8dexp 



r<t 



K 2 a 2 
t + a 



, a > 0, 



(3.84) 



and 



P{sup h* (s + r, u) — h* (s + t,u) > a} 

r<t 



< K 3 s~ K + K 4 (s + t) d exp 



t + a 



a>0. 



(3.85) 
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To prove (3.83) take a > 2\fdC\s. Since (x,u) < \\x\\2 < Vd\\x\\, as well as 
W^W < 1Mb < V^II^H, the left hand side of (3.83) is then bounded by 

2P{in V h (u, — C5«(s)) there is a B-particle outside C(a/Vd) 

at some time during [0, s] C [0, [2v^Ci] _1 a] } 

< 2P{nearest site a; to in S(u, —Csk(s)) with Na(x, 0— ) > 

lies outside C(a/(2y/d))} 

+ ^2 p. A^E or {number of S-particles outside C(a/(2\Q)) 

xeC(a/(2Vd)) 

at time [iVdCi^a] }, 

by an application of (3.7) and the argument following (3.25), very much as in (3.49). 
The first term in the right hand side here is at most 2exp[— K^a d ], and the second 
term is at most Kqot exp [-fiVdd]- 1 ^, by virtue of Theorem 1 in [KSb]. Thus 
(3.83) holds. 

The argument for (3.84) is basically already given in (3.19). Moreover, it is 
similar to, but simpler than the proof of (3.85) so we only prove the latter. If 
h*(s + t,u) = h, then all B-particles in V h (u, —C^k,(s + £)) are located in {x : 
(x,u) < h} at time s + t. If for some < r < t, h*(s + r, u) > h + a, then 
there is some S-particle p in V h (u, —Csk(s + r)) in {x : (x, u) > h + a} at time 
s + r. This p is also a particle present in V h (u, —C$k(s + t)) and even of type B 
in V h (u, — C$k(s + 1)) at time s + t, provided ||xo|| < CsK{s)/\fd (see (3.3)). Thus 
in this case p moved over a distance at least a/y/d during [s + r, s + t]. Therefore, 
the left hand side of (3.85) is at most 

P{\\x \\ >C 5 n(s)/Vd} 

+ P{some particle which starts outside C(3Ci(s + £)) becomes 

a S-particle in before time s + t} 

+ V ^^{sup \\S r -S t \\ > a/Vd} 
xec(3C 1 (s+t)) 

K q a 2 



<K 7 s- K + K 8 (s + t) d exp 



d(t + a) 



(see (3.48)-(3.50), as well as (2.42) in [KSa]). 

We can now proceed with subadditivity arguments. 
We introduce the random variables 

X(s) = [2VdC lS -h*(s,u)] + 

and the deterministic quantities Y(t) = 2\fdCiCQK{t), and let X'(t) be a copy of 
X(t) which is independent of X(s). Then (3.82) shows that, under (3.30), these 
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random variables satisfy 

P{X(s + t + C 6 K(t)) <(3}> P{X(s) + X'(t) + Y(t) <(3}- C s s~ K - 1 (3.86) 

for j3 > 0. Of course this inequality also holds trivially for (3 < 0. This is very close 
to the principal hypothesis of the lemma on p. 674 of [Ha] but we have to do some 
extra work because of the C^K(t) which appears in the argument on the left hand 
side of (3.86). From now on we take K = 4. Combining (3.86) with 

EX p (s + t + C 6 K(t)) 

^Vdd (s+t+C 6 fv(t)) 

= p / a p ~ 1 P{X(s + t + C 6 K(t)) > a}da 

Jo 

p2Vdd (s+t+C 6 K(t)) 

<p a p - 1 P{X(s) + X'(t)+Y(t)>a}da 

Jo 

+ C 8 [2VdC 1 ] p ( S + 1 + C Q K{t)) p s- K -\ 



we obtain 



EX(s + t + C 6 K(t)) 

< EX(s) + EX(t) + 2VdC 1 C 6 K(t) + K 10 [s + t + C 6 K(t)] S - K - 1 

(3.87) 



and 



EX 2 (2s + C 6 k(s)) < E[X(s) + X'(s) + Y(s)] 2 + K 10 [s + C^{s)] 2 S - K - 1 

< E[X(s)+X'(s) + Y(s)f + 4K 10 s- K+1 . (3.88) 

(3.88) holds for any s > sq, but so far (3.87) has only been proven under (3.30). But 
there is a simple replacement for this inequality that holds as soon as so < s < t. 
Indeed, assume that so < s < t, but tlogt > C-js 2 . We first observe that then 

X(s + t + C 6 K(t)) - X(t) < 2Vdd[s + t + C 6 K(t)] < bVdCtt, 

provided sq is large enough. Further, it follows from the simple inequality 

[a + b-c]+ -[a- d]+ < \b\ + [a- c]+ -[a- d]+ < \b\ + [c-d}~ (3.89) 

that 

X(s + t + C 6 K(t))-X(t) < 2VdC 1 [s + C 6 K(t)] + [h*(s + t + C 6 K(t),u)-h*(t,u)]-. 
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It follows that 

EX(s + t + C 6 K(t)) -EX(t) 

< 2\ / dC 1 [s + C 6 K(t)]+E\[h*(s + t + C 6 K(t),u) — h* (t, u)] ~ A h\fdC\t 

< 2Vdd[s + C 6 K(t)] + / P{h*(s + t + C 6 K(t),u) < h*(t,u) -a}da 

Jo 

< 2y/dd [s + C 6 K(t)] + bKaVddtt- 3 + 8d exp — ^— da 

J L s + C G K{t)+ai 

(by (3.84) with K taken as 3) 

< 2VdC 1 [s + C 6 n(t)} +K ll [s + C 6 K{t)] 1 / 2 < K 12 K(t). (3.90) 
Together with (3.87) this shows 

EX(s + t + C 6 K(t)) < EX(s) + EX(t) + K 13 k(s + t) (3.91) 
for all so < s < t. 

We shall next use a small variation on the argument of [Ha] to show that (3.91) 
implies 

X(u) := lim -Eh*(t,u)) exists. (3.92) 

t^oo t 

It suffices for (3.92) to show that 

lim -EX(t) = 2VdCt - X(u), (3.93) 

t^oo t 

because 

lim -E{h*(t, u);h*(t, u) > 2y/dC 1 t} = 0, 

t^oo t 

by virtue of (3.83). Now define for any M > e, 

to(M) = M, t k+1 (M) = 2t k {M) + C 6 K(t k {M)). 
Note that t k+ i/t k > 2, and hence t k {M) > 2 fc M, and for large k 



t k +i(M) C 6 /fclog3y/2 
2t k (M) ~ 2 V 2 k M 



) 1/2 , 



and for some K\±, independent of k > 0, 

fc-i 



1< yj t J±1 (M) _t k (M) K 14 

1 " 11 ~2^M) ~ ~M2F ~ 1 + M^- {3M) 



j 

Also, by (3.91), for all M > s + e, 

EX(t k (M)) < 2EX(t k - 1 (M)) + K 13K (t k (M)), k > 1 
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Consequently, 



EX{t k {M)) 
tk(M) 

EX(M) A 2^_i(M) ^K(t k -i(M)) -A- 2t J _i(M) 



<- 11 TTTaFS + K ^l^ 7, jm\ 11 

j=fc-£+i 



M l\ tj (M) f fc _,(M) V t,(M) 



< EX{M) [logM] 1 / 2 

In particular, for given £ > we can choose M > so + e so large, that 

ifis [log M] 1/2 M" 1/2 < e and EX(M)/M < liminf EX(s)/s + e. 

s^oo 

Then 

EXiUM)) ..EXis) 

\ fc ^. ;; < lim inf ^ + 2e, k> 0. 3.95 

Now let go > s o + be large. We shall expand qo as a sum of the form £*.(i) 
plus some error terms (see (3.97)) and obtain a corresponding bound for EX(qo) 
in (3.96). We define k(l) as the unique integer k for which t k < < £fc+i- We 
distinguish two cases. We are in the first case if go > + C6«(£fc(i)) + s o + -W- 
In this case we set qt = q - t k{1 ) - Ceftfei)) < -C6«(*fc(i)) = t fc (i)- 

Then s + M < qi < and 

EX(q ) < EX(t k{1) ) + EX( qi ) + K 13K (q ), 

by virtue of (3.91). If tfc(i) < qo < tk{i) + ^6^(^(1)) + s o + M, then, as in 
(3.89), (3.90), 

EX(q ) 

l>2\/~dC 1 qo 

< EX(t k{1) ) +2VdC 1 [q -t Hl) ]+ / P{h*(q ,u) -h*(t k(1) ,u) < -a}da 

Jo 

<EX(t k(1) )+K 16 K(q ) (by (3.84)) 

for a suitable large constant Kiq. If we are in the first case, we repeat the above 
procedure with qo replaced by q±. That is, we find k(2) such that t k ^) < Qi < £fc(2)+i 
etc. We continue to determine k(i) and qi until for the first time qi is in the second 
case, i.e., tfc(i+i) < Qi < *fc(i+i) + C6 K (tfc(i+i)) +so + M. Suppose this first happens 
at the index i - We then have 



EX(q ) < EX{t k{1) ) + EX( Ql ) + K 13 K(q ) 



< ... 



i + l io 

< Y,EX(t k(l) ) + (K 13 + K 16 )J2^)- (3-96) 

i=l i=0 
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Note that by construction, < t k ^ for 1 < % < io, and consequently, k(i+l) < k(i) 
for i < io- Therefore the above procedure ends at a finite io, and 

(K 13 + K 16 ) E <Q*) < (Kit + K ie) [ E K ^ + ^ Kl7 1-1° lo S «o] 1/2 ■ 

i=0 k:t k <q 

In addition we have either io = and qo > tfc(i) or zo > 1 and 

i io + l 

Qo =t k (i) + C 6 K(t k{1) )+q 1 = ■■■ = ^ [*fc(i) + C6«(*fc(i))] + Qi > E **(*)• ( 3 - 97 ) 

i=l 1=1 

Finally, we note that by definition of io, qi-i > sq + M, and therefore > M, 
for i < i . (3.96) and (3.95) now show that 

EX{qo) ^ 1^\ l^^EX(s) 



< - E c *« 



9o Qo U 



lim inf h 2e 



+ #17 

whence 



logQo] 1 / 2 ,.^ max j<M EX(j) 
+ l[tk(i +i) < M\ — 



Qo 



Qo 



EX(q) _ . r M(s) o 
lim sup < lim inf h 6e. 

q^oo q s ^°° S 

Thus the limit in (3.93) exists and we can use (3.93) to define \{u). 

We next turn our attention to the second moments. By definition < X(s) < 
2\fdC\s. With this inequality as a replacement of (11.10) and (11.12) in [Ha], we 
can start from (3.88) and imitate without essential changes the computations on p. 
676 of [Ha] or pp. 21, 22 of [SW] to obtain for any M > some 



oo 2 



E 



-l.V(/ A .(.\/)) 



(M2 k ) 2 



< oo. 



Since t k {M)/(M2 k ) > 1 (see (3.94)) we even have 

'jX(t k (l 
MM)} 



k=0 



and hence for any e > 0, 

OO ^ 
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By (3.83) also 

oo 

Y,P{X{t k {M)) ± 2^fdC 1 t k {M)-h*{t k {M) 1 u)} < oo, 

k=0 

so that for each fixed M > S4 and u £ S^ -1 , 
00 ^ 

J2 P {j^M)\ h *M M )^ u )- X ^)\ >e} <oo. (3.99) 

Of course (3.99) implies h*(tk(M),u)/tk(M) — > A(tt), almost surely. Since 
X(s) > by definition, 2^01 - X(u) > in (3.93), and hence X(u) < 2y/dd, 
as claimed. Finally, X(u) > C4 follows from Lemma 2 and the almost sure con- 
vergence of h*(tk(M),u)/tk(M) to X(u). In fact, (3.16) shows that almost surely 
h*(t k (M),u) = h(t k (M),u, -C 5 K(t k (M)) > C 4 t k (M) for all large k. 

Now choose a large Md > S4 and for some large integer r take Mj = Mo2 l / r , i = 
0,1,..., r - 1. Further take M r = *i(M ) and note that M i+1 /Mi -> 2 1 / 7 " as 
M — > 00 for fixed r and < z < r — 1. For given 77 > we can therefore first 
choose r large, such that 1 < 2 3 / r < 1 + 77, and then M so large that 

2 l/(2r) < Mj±i < 22 / r < i < r — 1. 
- Mi - - - 

By (3.94) we may further take Mq so large that 

2~^ — < **M < 2 V^, for M > M' > Mo, fc > 0. 
M> ~ t k (M') ~ M'' - - u ' - 

Once these choices have been made we take for {rij}j>o the collection of all dis- 
tinct £/o(Mj),0 < i < r — > 0, arranged in increasing order. Note that i 
only runs to r — 1 here. We claim that the collection {rij} in increasing order is 
{M ,Mi,...,M r _i,ti(Mo),...,ti(M r _i),rj 2 (Mo),...}. To verify this we merely 
need to check that tk(Mo) > tk-i(M r -i), since the other orderings are obvious 
from the monotonicity of tj(-). However, tk(M ) = tk-i{t\{Mo)) > tk-i(M r -i) is 
also easy from ^(Mq) > 2M > M r _i. This proves our claim. 
By construction we now have for all j > 0, 

2 i/(4r) < 2 -i/(4,) inf j t k {M i+1 ) fc>00<i<1 | < n i± i 

I t k (Mi) ~ ~ J - rij 

< 2 1 / r sup{ tfc( ^+ l) : k > 0, < % < r- l) < 2 3 / r < 1 + 77. 

1 tk ^ Mi > J (3.100) 

The leftmost inequality here shows that nj increases exponentially with j. 

It is simple to see that one can choose {rij(r}')} such that it contains the {n,j(rj)} 
from some index on, if rj' < 77, as claimed at the end of the Corollary. In fact if the 
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717(77) are constructed by the above method for some M ,r, then one can use the 
same construction for the rij(r)') based on Mg, r' with M' = t k (Mo) for some k and 
r' some integer multiple of r. 

Next, (3.81) holds, because by (3.99) 

00 

STp{\—h*(n k ,u)-\(u) >e\ 
k=o U Uk 

r— 1 00 ^ 

= ^^p{|_^/ i *(t fe (M i ),«)-A(7x)| >e} <oo. (3.101) 



i=0 k=0 

Thus also 

lim — h*(n k ,u) = X(u) a.s. (3.102) 

Now let < e < C 4 /2 < X(u) and ^4dC x < e/2. If 

— h*(n k ,u), ^—h*(n k+1 ,u)<2X(u)<4:VdC 1 (3.103) 
Tifc n fc+ i 

and 

h*(n k ,u) - K 18 K(n k ) < h*(t,u) < h*(n k+1 ,u) + K 18 K(n k+ i) (3.104) 

for all n k <t < n k +i, then, for k large enough and all n k <t < n k +i 

1 u*( \ ^ 1 u*t \ K is^{n k ) . 1 , 

— h (n k , u) — e < h (n k , u) < -h (t, u) 

n k n k+1 n k+1 t 

< — ft (n fc+ i,u) H < ft (n fc+ i,u) + e. 

n k n k n k+l ^ 105) 

By choosing K 18 large enough, and applying (3.84), (3.85), we can make 

00 

^^P{(3.105) fails for some t £ [n k , n k +i]} 

00 00 

< ^P{(3.103) fails} + ^P{(3.104) fails} 

k=0 fc=0 

00 00 

< 2 P { I — h * ( n k, - Mu) I > e} + 2 K * [n k ] " 4 

fc=o nfc fc=0 

00 

+ [ 8d + K 4 (n k+1 ) d j expi-K^K^r]- 1 \ogn k ] 

k=0 

< 00. (3.106) 

Since e > can be taken arbitrarily small, this, together with (3.102), proves the 
almost sure convergence in (3.80). The L p convergence along all reals in (3.80) 
follows from the almost sure convergence and the tail estimate (3.83). ■ 
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4. From half-space to full-space processes. 

The goal for this section is to prove that the B-particles in the full-space process 
do not spread faster than in the appropriate half-space process (see Corollary 8 for 
a precise statement). The first lemma establishes that for every u G S d ~ x there 
are deterministic vectors 14 such that for all r\ > there is with a probability 
close to 1 a S-particle in V h (it, — C$k{{1 + rj)rik)) "near" 14 at time njt, for all 
large k. Here rik is the rik(r]) of Corollary 5 and (14, it) has to grow essentially 
like h*(rik,u) ~ nkX(u) (see (4.1)). Apart from this growth condition the behavior 
of 14 as a function of k, u is unimportant for us. The only important aspect is 
that it is non-random, so that we can find, with high probability, a B-particle 
in a non-random location at which h*(n,k,u) is (almost) achieved. This will be 
used in the second lemma to concatenate V (u, —C^K^nk)) with another process 
which runs from time (1 + rj)nk to (1 + 77)71/5 + r k with r k also of order n k . By 
starting the second process at (14, (1 + v) n k) we will be able to assure that a 
B-particle at time (1 + rj)nk + r k in the second process is also a S-particle in 
V h (u, -C 5 K{{l + ri)n k + r k )). 

Lemma 6. Let u £ >S d_1 be fixed, and let n k = nk{rj) be as in Corollary 5. Then, 
for all < 77 < 64/(861) there exists a deterministic sequence of vectors {14} = 
{14(?7, «)} such that 

(V k (ri,u),u) = n k (rj)\(u), (4.1) 

and such that 

00 

^^P{in V h (u, —C^K{nk)) there is at time (1 + 77)71/- either no particle 

k=0 

at all in 14 + C(C2nkT]/4) or there is an A-particle in 14 + C(C2rikr]/4)} 
< 00. (4.2) 

Proof. Fix u E S^ 1 and e > 0. Let cr be a time which is so large that a > so (with 
so is as in Proposition 3) and such that 

-Eh*(a,u) - X(u) < (4.3) 

(see (3.27) for h*). Define the further times 

(7i = a, a j+1 = a + o-j + C 6 K(o-j), j > 1. 

Now apply (3.31) with the following choices: s = cr, t = aj, 7 = 7 S = Em*(a,u) 
(see (3.29) for m*) and 7 t = jEm*(cr,u). This yields 

P{G(a,0, (j + l)Em*(a,u),V h (u,-C 5K (a J+1 )),a J+1 ) 



> 



/ / P{h*(a,u) E dh,m*(a,u) E 7 + dm} 

J0<h<oo JmeR d 

P{Q(a - W - d,n - m,V h (u, -Csnfa)), a,)} - C 8 a 



K-l 
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provided (3.30) holds, that is, provided (aj + 1) log(oj + 1) < CVcr 2 . We start with 
j = r — 1, then use the case j = r — 2 with a, (3 replaced by a — h and (3 — <i, 
respectively, etc., all the way down to j = 1. With (/i*,m*), j > 1, i.i.d. copies of 
(/i*(cr, -u), m*(cr, w)) we obtain 

P{£ (a, /?, r^m* (a, u) , P ft (u, -C 5 K(a r )) , a r ) } 
> / / j]P{/i*(a, it) G dhj,m*(a, u) E 7 + dm.,} 

, ^ 7=1 



l<j<r-l l<j<r-l 

r— 1 r— 1 



x P\g{a -J2h 3 ,P-(r- l)d, 1 -J2 m Ji ph («» -C 5 k(<t)), a) } 
j'=i i=i 

- (r - ljCV"*" 1 
= pjin 7> h (u, -C 5 k(o-)) there is at time a a S-particle at some x with 

r— 1 r— 1 

(a;, it) + J^/i* > a and jjx^ + J^m* - 1 1 <&-{?- 
i=i j=i 

- (r - l)^"*" 1 

r r 

>P{J2 h *i * II EK - 7)|| < /3 - (r - l)d} - (r - l)^"*" 1 , 

i=i j=i (4.4) 

provided 

(<r r + 1) log(a r + 1) < C 7 a 2 . (4.5) 

It is easy to see by induction that each <jj is a continuous, increasing function 
of a on [0, oo). We further see by induction that > ka and <jj increases with j. 
Finally, we can for any fixed a > 1 find a K\ = K\{a) such that 

aK 1 2k{\ogk + 1) > C 6 k (a Kik 2 {log k + 1)), fc > 1, 

and <j\ < cri^i log 2. One more induction argument then shows that for all k > 1, 
< o--ftTi/c 2 (log fc + 1). Now let s > so be large and take r = [s 1 / 3 J . The preceding 
argument shows that we can fix a such that a r = s. Thus for j — 1 < r we have 
< cr r = s and cr^ < cr + (Tj_i +C6^(s). Consequently, s = cr r < ra + tCqk{s) = 
ra + Ls 1 ^ 3 JC'6k(s) = ra + o(s), and necessarily a ~ s/r ~ s 2 / 3 for large s. (4.5) is 
therefore automatically satisfied. If we further take 

a = ra[X(u) - -e], 

then, by (4.3) and the fact that Variance (h*) < K 2 a 2 (by (3.83)), 

P{ £ h* <a}<P{J2 W - ^J) < -rae/4} < |§. (4.6) 
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Further, fix S5 so large that 2se > roe > (1/2) se > 2rd ~ 2s 1 / 3 d for s > s 5 . Then 
we have similarly to (4.6), for s > S5 and (3 = se 

P{ II E K* - 7) || > /* - (r - 1)4 < P{ (I £ (mj - 7) || > rae/4} < ^. (4.7) 

The last two inequalities provide us with a lower bound for the right hand side of 
(4.4). We conclude that for s > s 5 

P{G (a, (3, rEm* (a, u) , V h (u, -C 5 kM) , a r ) } 

>l-^^-(r-l)C 8 .-->l--^ (4.8) 
(use any K > 1). Let 7^(77) be as in Corollary 5, and take s = n k = n k (t]). In 

1 /3 

agreement with our previous choice for r, a we then take r = \n k (rj)\ and a such 
that a r = n k (rj). Then, by going over to the complementary event in (4.8), we find 
for any rj > 0, that 

^P{in V h (u, -C 5 K(n k )) there is at time n k no S-particle 



k=0 



in F(n k [\(u) - ^-e],n k e,rEm*(a,u))} 



fc=0 n k e 

(recall that the rij grow exponentially). But (3.81) says in particular that 

00 

^P{in V h (u, -C 5 K(n k )) there is at time n k a B-particle 

k=0 

in r(n k [\(u) + -s],n k e,rEm*(a,u))} 

< 00. (4.10) 
We now take 

V k = V k (rj,u) =n k (r))\(u)u + rEm*(a,u). (4.11) 

Since m* is orthogonal to u (by definition), this choice of V k satisfies (4.1). More- 
over, (4.9) and (4.10) together give 



P{™ V h (u, -C 5K (n k )) there is at time n k no B-particle 

fc=0 

at any site x G V k + C(2\Qn k e)} 

00 

< £P{in V h (u, —C5K(n k )) there is at time n k no S-particle at any site x 



k=0 



with (x,u) e njfc[A(u) — -],nfc[A(u) + -] , Har 1 — rEm*(a, u)\\ < n k e) 
< 00. (4.12) 
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The convergence of the sums in (4.12) shows that almost surely, for all large 
rife (77), there is in V h (u, —C§K(nk)) a S-particle in V k + C(2\^dn k e) at time n k {r}). 
We claim that this implies that if we take e = 6*2/7/ (16c?) , then, in V h (u, — Cs^n*)) 
at time (1 + 77)71*;, all occupied sites in V* + C(C 2 nkf]/4:) are occupied by S-particles 
(and there are such occupied sites). More precisely, we claim that (4.2) holds. To see 
this we shall apply Lemma 4 with the following choices: s = n k , s = (1 + 77)77-/5, t = 
(1 + 77)77*; and finally y(n k ) is the location of any S-particle in V h (u : — 65^(77*;)) at 

time n* in the set Vk + C(C 2 rikr]/(8\/d)) , if such a S-particle exists. If several such 
.B-particles exist we pick the location of one of them according to some deterministic 
rule chosen in advance. On the event that no such S-particle exists we cannot apply 
Lemma 4, but this does not cause any problems, because (4.12) already tells us that 

00 

^P{no choice for y(nk) exists} < 00. (4-13) 

k=0 

If 7/(77,*;) exists, then there is automatically a particle in V h (u, —C^K,{rik)) at time 
77-*; at 7/(77-*;) 6 14 + C (C 2 n k i] / (8^/d)) . If this particle does not move a distance > 
C 2 n k 77/8 during [n*, (1+77)77*;], then it is in y(n k ) + C(C 2 n k r]/8) C V k + C(C 2 n k rj/4) 
at time (1 + 77)77*;. We recall further that all particles in V (u, — Cs^n*)) are also 
particles in V* . It follows that the k-th summand in (4.2) is bounded by the k-th 
summand in (4.13) plus 

P{||S„ fc „|| > C 2 n k r]/8} + P{B h (y(n k ), n k ;u, -C 5 «(n fc )) n K(y(n k ))} (4.14) 

(see (3.73) for JC(y)). The first probability in (4.14) is at most Kq exp[— ^777*77] by 
(2.42) in [KSa]. The last probability in (4.14) is by Lemma 4 at most 

P{y(n k ) i C(2C 1 77 fe )} + P{ (7/(77*), 77) < lc 4 77 fc } + n- k K ~\ (4.15) 

The first probability in (4.15) is 0[n^ K ~ 1 ^j by the estimates used for (3.59). 
The second probability in (4.15) is zero, because, by construction, y(n k ) e V k + 
C (C 2 n k r]/(8Vd)), so that 

(7/(77*0,77) > (V k ,u) -C 2 n k i]/8 

= UkX(u) — C 2 n k r]/8 > nk{C^ — C 2 i]/8) (see Corollary 5) > -C^n k . 

It follows that the sum of (4.15) over k is also finite, and this proves (4.2). ■ 

We can now show how to concatenate two processes as outlined before the last 
lemma. 

Lemma 7. Define 

H(t, 77) = h(t, 77, —00) 

= max{(x, 77) : x is occupied by a B -particle in at time t}. 

(4.16) 
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Assume that for some fixed u E 5' d_1 and \x > 

P{limsup -H(t, u)>n}> 0. (4.17) 

t^oo t 

Then 

AO) > fj,. (4.18) 

Proof. We divide the proof into 4 steps. Without loss of generality we assume 
H > 0. 

Step 1. For each small r\ > we choose 

ifi > > X(u) and if 2 = \= — • (4.19) 

ACWdKx 

We then define 

m k = m k {rj) = K 2 n k (r]), (4.20) 
where n k = n k {ri) is as in Corollary 5. We take 7] = i]o(e) > so small that 

li- ell 
l + Vo< ' 



Note that these definitions imply that for 1] < rjo, 

^±i = ^±i<l + r ? < »- £/2 . (4.21) 
m k n k ' t 1 - 3£/4 

Further, for small e > 0, define the events 
C k (n,fi-e) 

= {in there is a S-particle in the half-space S(u, m k (p — e)) at time m k } 
= {H(m k , u) >m k (fj,- e)}. 

In this step we shall show that for fixed e > and all < n < rjo(s), 

oo 

J2P{£k(v^-e)} = oc. (4.22) 

To prove this we shall show that 

P{C k {r], ii — e) occurs for infinitely many k} > 0. (4.23) 

(4.22) then follows from the Borel-Cantelli lemma. Now, (4.17) says that for every 
e > 

P{for infinitely many /c, H(t, u) > (fi — e/2)t for some t G [m k , m k+ i]} > 0. 

(4.24) 
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However, by (3.85) with h* replaced by H (this amounts to taking C5 = 00, which 
does not influence the estimate (3.85)) and with a = (e/4)mfc+i < (/i — e/2)mk — 
(ji - e)m k+1 (see (4.21)) 

P{H(t, u) > (fi — e/2)t for some t G [m k , m k+ i] but H(m k+ i, u) < (n — e)m k+ i} 
< P{ sup [H(r,u) - H(m k+ i,u)] > (/i - e/2)m k - (/i - e)m k+1 ) 

re[m k ,m k + 1 ] 

<K 3 (e,7 ] )[m k+1 ]- K . 

In particular, by Borel-Cantelli, the event in the left hand side here occurs almost 
surely only finitely often. Together with (4.24) this shows that 

P{for infinitely many /c, H(m k +i, u) > (fi — e)m k +i} > 0. 

This is the required (4.23). 

Step 2. The remaining steps are based on (4.22) only; (4.17) itself is not needed. 
With V k = V k (i], u) as in (4.11) we define an auxiliary process Q k = Q k {r], u) which 
is more or less the full-space process started at the deterministic space time point 
(Vfc, (1 + i])n k ). The only difference is that Q k only uses the particles which are at 
time in the "slab" 

{x : -n k /Ki < (x, u) - n k \(u) < Kin k }, (4.25) 

with Ki given by (4.19). Thus Q k is defined only from time (1 + rf)n k on. At 
time (1 + rf)n k it has at any x only the particles which started at time in the set 
(4.25). If no such particles exist, then there never are any particles in the process 
Q k . Otherwise, let z k be the nearest site to V k which is occupied at time (1 + r])n k 
by some particle, which at time was in (4.25). The types of all particles in Q k at 
time (1 + r])n k are reset to type A, except for the particles at z k , which are reset 
to type B. From time (1 + rf)n k the process then develops by our standard rules. 
Even though the process Q k is defined for all times in [(1 + rf)n k , 00) we are only 
interested in what happens during [(l + r])n k , (l + rj)n k +m k ]. Specifically, we define 
the events 

M k = M. k (r], \i — e) = {in Q k there is a S-particle in the half-space 

S(u, n k X(u) + m k (fi — e)) at time (1 + i])n k + m k }. (4.26) 
In this step we shall prove that 

J2P{M k } = oo. (4.27) 

To this end let us shift the event C k by (1 + rj)n k in time and by V k in space. Then 
C k goes over into the event 

C' k :— {in the full-space process started at (V k , (1 + rf)n k ) there 
is a .B-particle in the half-space S (-u, n k X(u) + m k (fi — e)) 
at time (1 + Vj)n k + m k }. 
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C' k \Mk can occur only if one of the following two events occurs: 

{at time (1 + rj)n k , some particle at the nearest occupied site to 
V k in the full-space process started at time outside the set (4.25)}, 

(4.28) 

or 

{in the full-space process started at (V k , (1 + v) n k) there is a 
particle which starts at time outside the set (4.25) and 
which coincides with a B-particle during [(1 + 77)77./-, (1 + rf)n k + m k ]} 

(4.29) 

(compare the argument for (3.47)). It follows that 

P{C' k \ M k } < P{(4.28) or (4.29) occurs}. 

But 

P{(4.28) occurs} 

< /-"{nearest occupied site to V k in V* at time (1 + rj)n k has 
distance more than log k from Vk } 

+ P{some particle which starts at time outside the set (4.25) 

is in Vk + C(K 4 log k) at time (1 + r))n k }. (4.30) 

Also, 

P{(4.29) occurs} 

< P{in the full-space process started at (V/-, (1 + rfjUk) there are S-particles 

outside Vk + C(2Cimk) at some time during [(1 + 77)77,*;, (1 + t])nk + rrik}} 

+ P{some particle which starts at time outside the set (4.25) 

visits V k + C{2C x m k ) during [0, (1 + v )n k + m k }}. (4.31) 

The first probability in the right hand side of (4.30) ) can be made 0[k~ K ) for 
any given K, by choosing K 4 large (compare (3.26)). The second probability in 
the right hand side of (4.30) is for large k no more than the second probability in 
the right hand side of (4.31). To estimate the latter, we merely point out that a 
particle which starts at some z outside the set (4.25) and visits V k + C(2Cim k ) 
during [0, (1 + r\)n k + m k ] has to move over a distance of at least 

||z-Vfc|| - 2C imk > d" 1/2 K(z - V k ),u)\ - 2C x m k 

= d~ 1/2 \(z,u) -n k X(u)\ -2dm k > n fc /(v / rfKi) - 2C 1 m k > n k /{2\fdK 1 ), 
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by virtue of our choice of m k . We leave it to the reader to use this to check that the 
last probability in (4.31) is 0{\n k ]~ K ) (see also the estimates in (3.24) and (3.50) 
or (3.76)). Finally, the first probability in the right hand side of (4.31) equals 

P{in V* there are S-particles outside C(2Cim k ) during [0, m k ]}, 

and this is 0([m k ]~ K ~ 2 ) , as in (3.25) and the lines following it. It follows from 
these estimates that P{C' k \ A4 k } < °o- In view of (4.22) and the fact that 
P{C' k } = P{£ k }, this implies (4.27). 

Step 3. In this step we show that 

P{Mk occurs for infinitely many k} = 1. (4.32) 

This is an easy application of Borel-Cantelli, because Ai k and Aig depend on 
particles which start at disjoint sets of sites (and are therefore independent) as 
soon as the set (4.25) and the corresponding set with k replaced by I are disjoint. 
If i > k, this is the case if n k {X{u) + Ki) < nt(X{u) — 1/Ki) and similarly if 
k > t. In particular, there is some integer K§ = K$(r]) such that M. k and A4g are 
independent as soon as \k — £\ > K^. Moreover, by (4.27), there is some integer 
j e [0, K 5 - 1] such that 

£ P{M k } = oo. 

k=j (mod K$) 

Thus (4.32) is true. 

Step 4. We now complete the proof of the lemma by showing that, almost surely, 
for all large k for which M. k occurs, also 

{in V h (u, — C 5 ((l + r])nk + rrik)) there is a S-particle in the half-space 

S(u, nkX(u) + mk{[i - e)) at time (1 + t])nk + rrik} 
= {h* ((1 + rj)n k + m k , u) > n k \{u) + m k {p> - e)} (4.33) 

occurs. This will indeed complete the proof, since we already know from Corollary 
5 that ((1 + i])n k + m fc )~V((l + rj)n k + m k ,u) -> X(u). Thus (4.32) and (4.33) 
will imply, for all e > 0, < r] < rjo(e), 

X(u) > liminf ^ X(u) + — ^ (u - e) 

fc^oo L(l +7])n k + m k {l + i])n k + m k 



1 + T] + K 2 l+T] + K 2 

and hence 



X(u)>^-(p-e). (4.34) 
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Now to prove (4.33), we write, as in the lines following (4.25), Zk for the nearest 
site to Vfc at time (1 + rf)nk which is occupied by a particle which started at time 
in (4.25). We already proved that, almost surely, (4.28) occurs only finitely often. 
Thus, except for finitely many k, Zk actually equals the nearest occpied site to at 
time (l+r])nk inV^. Since the set (4.25) is contained in S(u, 0) C <S(w, — C^K{nk)), 
Zk is also the nearest occupied site to "14 at time (1 + rf)nk in V h (u, —C^K{nk))- 
By virtue of Lemma 6, we further know that, a.s. for all large k, Zk is occupied 
by .B-particles at time (1 + r\)nk in V h (u, — Cs(nfc)) for all large k. By using 
the monotonicity property of Lemma C we conclude that, almost surely, for all 
large k all the B-particles in Q k at time (1 + rj)nk + mk are also B-particles in 
V h {u, —C§k{{1 + rj)rik +mfc)). In particular, 

h* ((1 + r])n k + mk, u) > n k X(u) +m k (p- e) 

for all large k for which Mk occurs. This is the required (4.33). ■ 
Corollary 8. For every unit vector u 

lim -H(t, u) = X(u) almost surely and in LP for all p > 0. (4.35) 

t^oo t 

(t runs through the reals here). Moreover, for nk = nk{rf) as in Corollary 5, we 
have for any 5 > and rj > 0, 

X . 

J2 P i\— H ( n k,u) - X(u)\ > 6} < oo. (4.36) 

k=o Uk 

Proof. By the monotonicity property of Lemma C 

H(t,u) > h*(t,u) on the event {||x || < C 5 K(t)/Vd} (4.37) 
(see also the lines after (3.48)). Thus, by the estimate (3.26) 

liminf -H(t, u) > lim -h*(t,u) = X(u) 

t^oo t t^oo t 

(see Corollary 5). In the other direction, we have from Lemma 7 that 

1 

Pjlimsup —H(t, u) > fx} = for all fx > X(u). 

i^oo t 

This proves the almost sure convergence in (4.35). The LP convergence follows from 
the almost sure convergence and the tail estimate 

P{H(s, u)>a}< exp[-Kiai\ for a > ^Vdds, s > s 3 , (4.38) 
which can be proven in the same way as (3.83) (or we can take C5 = 00 in (3.83)). 
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As for (4.36), we have by (4.37), (3.81) and an estimate like (3.26) that 

J2 P {— H ( n k,u) < X(u) -5} < oo. (4.39) 



fc=0 



For the other direction, we begin with an indirect argument. Assume, to derive a 
contradiction, that for some 5 > and < r\ < C^K&Ci) 



- 1 

J2 P i — H ( m k,u) > X(u) + 5/2} = oo, 



with m k = m k (rj) as in (4.20). This is just (4.22) with \x — e replaced by X(u)+S/2. 
By steps 2-4 of the proof of Lemma 7 we then have that (4.34), again with [i — e 
replaced by X{u) + 5/2, holds. This is impossible for r\ < K 2 5 / (2X{u)). Thus for 
all 5 > 0, < r] < C 4 /(8Ci) A K 2 S/(2X(u)), it is the case that 

oo 

V?{ H(m k ,u) > X(u) +6/2} < oo. (4.40) 

k=o K 

Finally, for given k, let i = £(k) be determined by mi < n k < me+i- We now use 
that 

P{H(n k , u) > n k (X(u) + 5) } < P{H(m e+1 , u) > m e+1 (X(u) + 5/2) } 

+ P{H(m e +u u) - H(n k , u) < m e+1 (X(u) + 5/2) - n k (X(u) + 5)}. 

(4.41) 

But, by (3.84) (with C 5 taken to be infinity) we have 

P{inf H(s + r,u)- H(s, u) < -a} < K 3 s~ K + 8dexp [ - ? a > o. (4.42) 

r<t L t + a J 

Moreover, m^+i < (1 + r])mi < (1 + rf)n k (see (4.21)). Therefore the second term 
in the right hand side of (4.41) is at most 

P{H(m e+1 , u) - H{n k , u) < n k [(l + rj) (X(u) + 5/2) - (X(u) + 5)] < -n k 5/4} 
< K z n- k K + K 6 e W [-K 7 n k 5 2 /( V + 5)}, 

provided 



C A K 2 5 

2Xju) ' 4(A(u) + 5/2) 
It follows that under this last condition 



71 < mm W' 2XUTV 4(X(u) + 5/2) \- (443) 



J2P{H(n k ,u) >n k (X(u) + 5)} 

k=0 

oo 

< P{H(m m+1 ,u) > m m+1 (X(u) + 5/2)} + 0(1). 



k=0 
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The right hand side here is finite by virtue of (4.40), because m^ k ) = K 2 n^ k ) < 
n k < K2ri£(k)+i forces \£(k) — k\ < Kg for some Kg which is independent of k (see 
(3.100)). Finally, we may drop the condition (4.43), because if r\ does not satisfy 
this condition, but rj' does satisfy this condition, then we may choose {n k {j]')} so 
that it contains the tail of {n k {r])}, by Corollary 5. By this inclusion and by what 
we just proved 

oo 

£ P{H(n k (rj), u) > n k (rj)(X(u) + 5)} 

k=0 

oo 

< £ P{H(n k (rf),u) > n k ( V ')(X(u) + 5)} + 0(1) < oo. ■ 

fc=0 

5. Proof of the shape theorem. Now that we have shown that the spread 
of the .B-particles in the full space process has a definite speed in each direction, 
the half-space processes are no longer of importance. In fact Corollary 8 contains 
Theorem 1 in the one- dimensional case (with Bq = [— A(ei), A(ei)]). For the higher 
dimensional case, we shall show in this section how to go from the existence of 
lim t _ ! . 00 (l/f).ff"(t, u) for all u G S^ -1 to the full shape theorem. This should work 
for a fairly general class of processes. The idea to derive the shape theorem via 
results on the propagation of half-spaces we learned from [GG]. However, the details 
in our case differ from those in [GG]. 

The remaining problem in dimension d > 1 is that even if we know that H(t,u) 
grows at rate X(u), it only tells us that there exist B-particles at time t at some 
random site xt for which (xt,u) ~ t\(u). It does not tell us where the points xt 
near the hyperplane {x : (x,u) = tX(u)} are. In particular, it does not guarantee 
that we can find Xt which converge in direction to a prescribed unit vector, i.e., for 
given v G S' d_1 we do not know whether we can choose Xt such that Xt/||xt||2 — > v. 

To attack this problem we first write down the conjectured limiting shape B 
in terms of the function A(-) on 5 >d_1 . This conjectured B is convex (for trivial 
reasons). We then show that we can guarantee a^t / 1| cc* || 2 — ► v if v corresponds to 
a so-called exposed point of the convex set Bq. Using some further properties of 
convex sets, as well as approximate convexity properties of the set of points which 
can be reached by the S-particles in a large time, we can then show that the limiting 
shape result (1.3) holds. 

The convergence result (4.35) suggests that the limit set B in (1.3) should be 
given by 

B = {z E R d : (z, u) < X(u) for all u E S^ 1 }. (5.1) 

Clearly this set Bq is a closed convex set. In fact it is also bounded and hence 
compact, because A (it) < 2\fdC\ for all u. The origin is an interior point of So 
because X(u) > C4. We call a point w G dBo an exposed point of Bq if there exists 
a supporting hyperplane {z G M d : (a, z) = b} of B which contains w, but no other 
point of B . Thus 



(a, w) = b but (a, z) < b for all z G Bq \ {w}. 



(5.2) 
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Note that this forces a^O. We now show that V* indeed grows in the direction 
of an exposed point at the rate which is necessary for (1.3). 

Lemma 9. Let w be an exposed point of B and let (a, b) £ W 1 x R satisfy (5.2). 
Let u = a/||a||2. Then, there exists a sequence e n | such that 

P{Af n (w, e n ) occurs for all large integers n} = 1, (5-3) 

where 

J\f n (w,e) := {in there are at time (l + 8e/C2)n occupied 

sites in nw + C(2en) and all these sites are in fact 
occupied by B-particles at time (l + 8e/C2)n}. (5.4) 

Also, define 

O n (w, 5) = [in there is at time n a B-particle in nw + C(5n)}. 
Finally, let = nk(t]) be as in Corollary 5. Then for all 5, n > 

oo 

^[l-PR lW K^)}] <oo. (5.5) 

k=0 

Proof. Order the vertices of Z d in some deterministic way, for instance in the lex- 
icographic way. Let xt be the first vertex x in this order which is occupied by a 
-B-particle in V* at time t and with (x,u) = H(t,u). By (4.35), almost surely, 

1 1 

-(xt,u) — >■ X(u) = lim -H(t,u) (5.6) 

t t^oo t 

as t — > oo. Moreover, by (4.36), for each S > 0, r] > 

CO 

J2P{\ — (x nk ,u) - X(u)\ >5}<oo. (5.7) 
k=o Uk 

We want to show that for each 5 > 

-Pj II— x n — w\\ < 5 for all large integers n\ = 1. (5.8) 
Note that w E B implies 

(w,u)<X(u). (5.9) 

Recall next that P{x n (£ C(2dn)} < K & n- K - d ~ 1 , by virtue of (3.49) or the 
estimates for (3.25). So, 

P{x n E C(2C 1 n) for all large n} = 1. (5.10) 
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Also 



J2 p i x n k £C(2Cm k )} < OO. 



(5.11) 



k=0 



So, we can ignore the events {x n £ C{2C\n)}. 

Next, let v G >S d_1 be a unit vector which is not a multiple of w. We claim that 
there exists some 5 = 5(v) > such that 



— v 



Fn 2 



< 5 i.o. 







and 



fc=0 



Fn fc ||2 



< < OO 



(5.12) 



(5.13) 



(i.o. stands for infinitely often). To prove this, note first that (5.12) holds if 
(v,u) = 0, because 



liminf(- 



Fn 2 



-,u) > 



X(u) 



> 



C 4 



limsup n ^ 00 \\x n h/n ~ 2Vdd 



by virtue of (5.6), (5.10) and the fact the X(u) G [C4, 2\/~dC\]. Similarly, (5.13) 
holds if (v,u) = 0, by virtue of (5.7) and (5.11). To take care of other vectors v, 
define for any y G R d \ {0}, with (y, u) ^ 0, 

y = the unique multiple of v which satisfies (y, u) = 6/||a||2- 

In particular, y lies in the in the supporting hyperplane {z : (a,z) = b}. Now, by 
assumption v 7^ w, so that v £ B . By definition of B this means that there exists 
some u' G S^ -1 such that (v,u') > X(u'). We can then find S > and i] > such 
that (z, u') > (1 + ri)X{u') for all z G with \z - v\ < 5. Thus, if 



Fn 2 



<5, 



then 

(Xn,u') = ((x n /\\x n \\ 2 ),u') > (l+7])X(u). 

In addition, by (5.6) and (5.9), 

lim —(x n ,u) = X(u) > (w,u) = - — — (see (5.2)), 
n^oo n a h 



(5.14) 



while, by definition of y, 
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Moreover, we must have 

||a|| 2 (w, u) = b > (5.15) 

by (5.2) and the fact that G B . Consequently, x n = 7 n :r n for some reals 7 n 
which satisfy 7 n /n — > 1. Together with (5.14) this would imply 

(x n ,u') > n(l + rj/2)X(u') 

for large n. But, P{(x n ,u') > n(l + rj/2)X(u') i.o.} = 0, by virtue of (4.35) with 
u replaced by u' and the fact that H(n,u') > (x n ,u') (by definition of H). Thus 
(5.12) holds for the chosen 5. Similarly, (5.13) follows by means of (5.7) with v! 
instead of u. 

Now, for any e > the compact set 

W(e) : = {ze S*- 1 : z = -f^- for some x G C(2Cm) with 



F 2 



w 

(x, u) > n\(u)/2, \\z — - — — 1| > e} 

\\ w \\2 



= {zeS d ~ l :z = tt^tt for some x G C(2C t ) with 
iFlb 

w 

(x,u) > \{u)/2, \\z- T r ^ r \\ > e} 

is independent of n and is covered by finitely many neighborhoods Ui, . . . , Un of 
the form U t = {z E : \\z - v t \\ < 8(vi)} with v t G S d -\ Thus, by (5.12), 

P{x n /\\x n \\ 2 G W(e) i.o.} = 0. This holds for all e > 0. In view of (5.6) and (5.10), 
this implies 

MFnlh H^lb J 

In turn, this together with (5.6) implies 

n\(u) (x n ,u) (w,u) 
hm j — = hm -j — = - — — a.s. 

n^oo ||x n || 2 n->oo \\x n \\2 \\w\\2 

Since (w, u) ^ (see (5.15)), ||x n ||2 ~ n||w||2A( , u) / (w, u) and 

hm — x n = a.s. (5-17) 

n^oo n u) 

To complete the proof of (5.8) we show that 

\{u) = {w,u). (5.18) 

Indeed, we already saw that {x n , u) = 6/||a||2 = (w, u). We also saw that x n = 'jnXn 
with 7 n ~ n. Therefore (x n /n,u) ~ (x n ,u) = (w,u). On the other hand, (5.17) 
implies that lim n ^ 00 (x n /n, u) = X(u). Thus (5.18) and (5.8) hold. 
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We now also obtain (5.5). Indeed, essentially the same argument as for (5.16), 
but now using (5.13) instead of (5.12) gives 

oo 

> s\<oo. (5.19) 

Consequently also 

E nn 1 , \ \\ x n k lb , \ i SVd u .. > 

, _ [ n k n k \\w 2 n fc 

fc=0 

Together with (5.7), (5.18) and (5.11) this last relation yields 

oo II II 

V PUXOu) - ,?* „ A(m)| > dddd] < oo. 
f-1 n k \\w\\ 2 

k=0 

Thus, for a suitable constant K? 

OO II || 

^pjl^J _ J > ^ < qo _ 

fc=o nfc 

Together with (5.19) this finally gives for some other constant Kg 

oo 

52[l-P{O nhM (w,K 8 S)}] 

k=0 

oo 



^E p {ll^r-HI>^}<~- 



Since this holds for any 6 > 0, this is equivalent to (5.5). 

The preceding (see (5.8)) shows that there exists a sequence e n — > 0, and random 
vertices x n such that with probability 1, for all large n, 

x n G nw + C(e n n) and (x n , n) occurs, (5.20) 

where 

t$f(x,s) := {there is B-particle at x at time s in ?^}. 

Now take 

n:=n(l + ^) 

and define the event 

lZ(x, n) = {there is some particle in which lies in x + C(C 2 (n — n)/2) 
= x + C(Ae n n) but is of type A at time n}. 



SHAPE THEOREM FOR SPREAD OF AN INFECTION 



55 



We shall complete the proof by proving that the event 

{for infinitely many n there exists an x n for which 

B^(x n , n) n 7Z(x n , n) occurs} (5-21) 

has probability 0. First we show that this will indeed prove the lemma. The 
probability that any particle which is in nw + C(e n n) at time n is outside nw + 
C{2e n n) at time n is bounded by 

K 9 [e n n] d P{ sup \\S r \\>e n n}. (5.22) 

r<n—n 

Without loss of generality we can let e n go to so slowly that for large n this 
expression is no more than n~ K ~ x (by (2.42) in [KSa]) and such that 

e n > n~ 1/2 . (5.23) 

From this and the fact that B^(x ni n) occurs for all large n, we conclude via the 
Borel-Cantelli lemma that almost surely, for all large n there are particles in in 
the set nw + C(2e n n) at time n. The fact that (5.21) has probability will then 
imply that 1Z(x n ,n) must fail for all large n. But this implies that a.s. there are 
particles in V* which lie in nw + C(2e n n) C x n +C(4e n n) at time n, and all of these 
particles must have type B. This is the desired result (5.3). 

It remains to prove (5.21). But this is almost immediate from Proposition 1. 
Indeed, 

P{B f (x n ,n)Ml(x n ,n)} 

< P{x n C(2Cm)} + ^2 P{& f ( x , n ) but there is a particle in V f of 

xEC(2C 1 n) 

type A at some z G x + C(4e n n) at time n} 

< KQn~ K ~ d ~ 1 + P{x is occupied at time n in and in the full-space 

xeC(2C\n) 

process started at (x, n) there is an A-particle 
at some z G x + C(4e n n) at time n}, 

where we used Lemma C for the last inequality. As in the estimate for /C2 in (3.77), 
by (3.9) and (2.4) with K replaced by 2K + 2d, the last sum here is at most 

K 10 n d (Mt hand side of (2.4) with t = n-n = 8e n n/C 2 ) < K xl n~ K 

(see (5.23)). ■ 

The preceding lemma shows that the set B(t) grows in the direction of the 
exposed points of Bq in 8Bq at the "right" speed. More specifically, if w is such a 
point, then almost surely, for all large t, there exist points w(t) G (l/t)B(t) such 
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that w(t) — > w. We merely have to choose n in Lemma 9 such that n(l+8e n /C2) < t 
but n/t — > 1, and then w(t) a point in B(n) D [nio +C(2e n n)]. Lemma 9 guarantees 
that this last intersection is nonempty for large n. The next two lemmas will 
show that the same is true for any point w G OBq. This is basically done by 
concatenating a number of paths which produce B-particles at c^mu^ for exposed 

points w n ,i and Yl,i=i a i w n,i w,oti > 0, X)i=i a i = 1- Lemma 10 contains the 
basic technical step. It explains how the concatenation works; this is basically the 
same construction as in the proof of Lemma 7. 

Lemma 10. Let u>i,W2 G dB . Assume that there exist e n > such that e n — > 
and such that (5.3) holds with w replaced by w\, that is, 

P{N n {wi,e n ) occurs for all large integers n} = 1. (5-24) 

( We are not assuming that w\ is an exposed point of B .) In addition, assume that 
for all 5, rj > 

CO 

J2^-P{O nk (v)( w 2,S)}} <oo (5.25) 

k=0 

(see Corollary 5 for nk = n,k(r})). Let < a < 1 and rj > 0. Then there exist 
5 n > such that 5 n — > and swc/i t/iat 

-P{-A/" n (awi + (1 — a)w2, <5 n ) occurs for all large n} = 1. (5.26) 



Proof. Fix < a < 1. Also fix 



5 > and < rj < 5/2 



for the time being. Take 



Pk =Pk(v) = 



a 



1 — a 



n k {ri) 



and 



Qk = Qk(v) = (1 + 8e Pk{7l) /C 2 )pk- 



Define O' (^2, 5) as the translate by (pk(v) w ii Qk(v)) ( m space-time) of O nk (w2, S). 
Explicitly, 

0' nk (w2,5) = {in the full-space process started at (pk(v) w ii Qk(v)) there 
is at time qk + nk a S-particle in PkWi + UkW2 + C(Sn} t )}- 

(We suppress the dependence on w\ and r\ in this notation). Also let 

Zk = nearest occupied site to PkWi in at time qk- 



SHAPE THEOREM FOR SPREAD OF AN INFECTION 



57 



Since P{0' nk (w2, 5)} = P{O nk (w2, 5)}, assumption (5.25) implies that almost 
surely, 

0' nk (w2,S) occurs for all large k. (5.27) 
Also, by assumption (5.24), almost surely, 

Af Pk (wi,E Pk ) occurs for all large k. (5.28) 

Now consider a k for which J\f Pk {w\ , e Pk ) fl (D' nk (w2,S) occurs. By the definition of 
J\f Pk this implies that z k lies in p k wi + C(2e Pk p k ) and that the particles at z k at time 
q k have type B in V? . Therefore the resetting of the types to start the full-space 
process at (p k w±, qk) does not change the type at z k . By the monotonicity property 
of Lemma C, V? therefore has at least as many S-particles at any space-time point 
(x,t) with t > qk as the full state process started at (p k wi,q k ). Since 0' nk (w2,S) 
occurs this implies that in V? there is a S-particle in PkWi + n k W2 +C(5rik) at time 
Qk + ™fc- 

Let the nearest B-particle to p k wi + n<kW2 in V? at time qk + rik be at the 
position y k , so that B^(yk,qk + n k ) occurs. The last paragraph gives us that 
\\yk —PkW\ —nkW2\\ < Snk- These are only statements for the times qk +nk- Since 
(5.24) requires that certain events happen for all large n we now first show how to 
go from the qk + rik to general integers n. For any large n let k(n) be such that 
qk + rik < n < qk+i + rajfe+i- Then for large n 

qk+rik <n< (q k + n k ){l + 2rj) < (q k + n k ){l + 5), 

since n k+ i/n k < 1 + r\. Also by our choice of p k , qk 

\\p k wi + n k W2 - n[awi + (1 - a)w 2 ] \\ < K 12 Sn. 

Thus, on (yk, qk+nk), there is a S-particle at y k G n[awi + (l — a)w2]+C((K 12 + 
l)Sn) at time q k +n k . Moreover, as in (5.22) we have 

P{in there is a S-particle in n\aw\ + (1 — a)w2] + C[(Ki 2 + l)Sn) 

at time q k + n k which is no longer in n[awi + (1 — a)w2] + C[(Ki 2 + 2)5n) 
at time n} 

= 0{n- K ). 

Thus, almost surely, there is in V? for all large n a S-particle in n[awi + (1 — 
a)w2] +C((K 12 + 2)5n) at time n. We can now proceed as in Lemma 9. Essentially 
as in (5.21) and in the lines following it we now have that almost surely 

{there is some y G n[aw\ + (1 - a)w 2 ] + C((K 12 + 2)Sn) for which 
B^(y,n) occurs, but in there are either no particles or an A-particle in 
n[a Wl + (1 - a)w 2 ] + C(2(K 12 + 2)8n) at time (l + 8(Ki2 + 2)5/C2)n} 

(5.29) 
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occurs only for finitely many n. This shows that 

P{N n (awi + (1 — a)w2i (K12 + 2)5) occurs for all large n} = 1. (5.30) 

This holds for all 5 > and rj < 5/2. However, (5.30) is already independent of rj, 
so that it holds for all 5 > 0. There then also exists a sequence 5 n — > such that 
M n {oLWi + (1 — ol)w2-> 5 n ) occurs almost surely for all large n. ■ 

Proof of Theorem 1. We shall prove (1.3) with the Bq defined in (5.1). For 
the right hand inclusion in (1.3) we note that for any e > there exists finitely 
many halfspaces {z G M. d : (z, Ui) < A(-Ui)}, 1 < i < N, with Ui G such that 

N 

f]{z G R d : (z, Ui ) < X( Ui )} C (1 + e/3)B . (5.31) 

Indeed, B is contained in the cube C := fl^-fz G M. d : — A(e^) < (z, e^) < A(ei)} 
(with ei = i-th coordinate vector), and by compactness, C\(interior of (1 + s/3)Bq 
is covered by finitely many relatively open subsets of C of the form C n {z G M d : 
(z, tt) > A(w)}. In addition to (5.31) we know from (4.35) that, almost surely, 
H{t, Ui) < t(l + e/3)\(ui) for all large t and i = 1, . . . , N. Consequently, almost 
surely 

N 

B(t) Ct(l + e/3)f){z eR d : (z, Ui ) < X(u t )} C (l+e/3) 2 tB 
i=i 

for all large t. Thus the right hand inclusion in (1.3) holds. 

For the left hand inclusion in (1.3) we first observe that by Lemma 9, the hy- 
potheses (5.24) and (5.25) of Lemma 10 hold for all exposed points 1^1,^2 G dB . 
It then follows from Lemma 10 that (5.26) holds. In turn, (5.26) states that the 
hypothesis (5.24) with w\ replaced by aw\ + (1 — a)u>2 is satisfied. Therefore, if 
103 G 8Bq is also an exposed point of So and < (3 < 1, then we get from Lemma 
10 that there exist 5' n — > such that 

P{Af n (f3awi + (3(1 — a)u>2 + (1 — f3)w3, 5' n ) occurs for all large n} = 1. 

But as a and (3 vary over (0, 1), (3aw\ + (3(1 — a)u>2 + (1 — (3)ws varies over the 
convex combinations ai?«i + 0:2^2 + &3W3 with ctj > 0, Yli=i a i = 1- ^ e can 
repeat this procedure to obtain that for each convex combination X^i=i a i w i with 
oti > 0, Ylfi=i a i = 1 an d Wi G 8Bq exposed points of So, there exist 5 n — > such 
that 

k 

P{J\f n ( ^ otiWi, 5 n ) for all large integers n} = 1. 

i=i 
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P{in V* there are at time (1 + 85 n /C2)n B-particles in 

k 

n otiWi + C(2nn) for all large integers n} = 1. 

i=i 

In turn, this means that if for a given vector v and rj > we can find ct^tUj as 
above such that \\v — X^=i Q^ill — Vi then also 

P{in there are at time (1 + 85 n /C2)n B-particles in 

nv + C(3rjn) for all large integers n} = 1. (5.32) 

If v is such that there exist < oo, a\ > and G 8Bq exposed points of Bq 

such that Yli=i a i = 1 an d — X)i=i a ^ w ^ II ~" ^ (as r — > oo), then (5.32) holds 
for each r\ > 0. For such d there then exist r\ n — > such that almost surely, for all 
large n there exist S-particles within distance 4?y n n of nv at time (1 + 8S n /C2)n, 
for some 5 n — > (S n and ?y n may depend on v). 

The last statement applies to each v & Bq, because each such v is a convex 
combination of at most (d + 1) extreme points of So (see [Ru], Theorem 3.22 
and Lemma following Theorem 3.25) and the exposed points of Bq are dense in 
the extreme points (Strascewicz' theorem; see Theorem 18.6 in [Ro]). Thus, by 
applying the last result to a fixed v & Bq with n = [(1 — e)t\ and < e < 1, we 
find that almost surely for all large t, 

at time (1 + 85 n /C 2 )n there exists a site v n with 

\\v n — nv\\ < 4i] n n, which is occupied in by S-particles. (5.33) 
We claim that 

P{(5.33) holds, but not all sites in (1 - e)tv + C(C 2 et/A) belong to B(t)} 
< K 13 t~ K . (5.34) 

This is an easy consequence of (3.9) and Theorem A. Indeed, from (3.9) with (X, s) 
taken to be (v n , (1 + 8d n /C2)nj and 

A = {not all vertices in C(C2en/2) have been visited by a S-particle by time en/2} 
= {C{C 2 en/2) <f_ B(en/2)}, 

we see that the probability in (5.34) is for large t at most 

K 14 n d P or {C(C 2 en/2) <£ B(en/2)} < K 15 n~ K < K 13 t~ K 

(for the first inequality here we used Theorem A with K + d in the place of K). 
This establishes the claim (5.34). 
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To obtain Theorem 1 we now choose for a given e a finite number of vectors 
. . . , v^ N ' in So such that each v & Bq satisfies \\v — v^ r ' \\ < C2E/4 for at least 
one r. This means that 

BoC U [v^ +C(C 2 e/4)]. 

l<r<N 

Moreover, by (5.33) and (5.34) it holds almost surely for all large t that 

|J [(1 - e)tv {r) + C(C 2 et/4)] C B(t). 

l<r<N 

Together, these last two inclusions imply that almost surely the left hand inclusion 
in (1.3) holds for all large t. 
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